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ABSTRACT 


1.  PROGTEST.  This  Computer  System  for  the  Analysis  of  Computa¬ 
tional  Computer  Programs  consists  of  three  interrelated  subsys¬ 
tems: 

POINTCOMP  consisting  of  the  major  components  POINTCOMP 
and  VARVARYl 

ANALYZ  containing  the  ANALYZ  routine 

GRID  consisting  of  the  GRID,  REARRANGE  and  DIFFQUOT 

routines  as  major  components 

2.  POINTCOMP 

a.  The  POINTCOMP  routine  of  the  POINTCOMP  subsystem  utilizes 
the  method  of  steepest  ascent/descent  to  generate  input  variable 
values  maximizing/minimizing  the  output  of  the  program  being 
tested  (one  output  at  a  time).  The  points  generated  by  the  rou¬ 
tine  are  printed  out  and  may  be  used  by  the  analyst  to  determine 
whether  increasing  or  decreasing  the  variable  increases/decreases 
the  output,  or  whether  a  maximum/minimum  is  reached.  The  routine 
also  prints  out  sensitivity  information  usable  by  the  analyst  in 
determining  the  comparative  effect  of  each  variable  upon  the  out¬ 
put  variable.  Points  generated  by  this  routine  are  used  as  input 
by  some  of  the  other  subsystems. 

b.  The  VARVARYl  component  of  the  POINTCOMP  subsystem  uses  in¬ 
put  variable  values  produced  by  POINTCOMP  or  chosen  by  the  ana¬ 
lyst.  For  each  combination  produced  or  chosen  by  the  analyst, 
each  variable  is  varied  uniformly  through  a  range  while  keeping 
the  other  variables  fixed.  The  tested  routine  is  evaluated  at 
each  new  point  and  the  gradient  is  computed  at  each  new  point.  In 
addition,  statistics  showing  the  average  marginal  return  over 
various  portions  of  the  variables  range  are  printed  out.  The  out¬ 
put  of  this  program  can  be  used  to  generate  curves  showing  the  ef¬ 
fect  of  the  variable  being  varied  upon  the  output. 

3.  ANALYZ.  The  ANALYZ  subsystem  uses  a  file  of  gradients  as  in¬ 
put.  This  type  of  file  is  produced  by  several  subsystems.  The 
subsystem  computes  and  outputs  sensitivity  statistics  derived  from 
the  gradients. 
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4.  GRID.  The  GRID  routine  of  the  GRID  subsystem  utilizes  a  user 
defined  grid  to  intensively  analyze  a  small  area  of  input  values. 
The  subsystem  evaluates  the  routine  to  be  tested  at  each  node  in 
the  grid,  and  prints  out  the  input  values,  the  evaluated  values, 
and  the  gradients  at  each  point.  The  REARRANGE  routine  of  the 
GRID  subsystem  rearranges  the  data  to  show  the  variable  by  vari¬ 
able  variation,  and  the  DIFFQUOT  routine  computes  difference  quo¬ 
tients  for  the  generated  points  for  one  variable  at  a  time. 


iv 
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PROGTEST:  A  COMPUTER  SYSTEM  FOR  THE  ANALYSIS  OF 
COMPUTATIONAL  COMPUTER  PROGRAMS 


CHAPTER  1 


OUTLINE  OF  SUBSYSTEMS 


1-1.  SUBSYSTEM  DESCRIPTIONS 
a.  PQINTCOMP  Subsystem 

(1)  PQINTCOMP.  The  PQINTCOMP  routine  of  the  POINTCOMP  sub¬ 
system  creates  input  variable  values  by  using  maximum  ascent/de¬ 
scent  techniques  to  maximize/minimize  the  output  of  the  program 
being  tested.  The  output  from  POINTCOMP  includes: 

(a)  A  file  in  standard  format  for  which  data  lines  are 
comprised  of  variable  values  and  the  corresponding  value  of  the 
program  being  tested.  This  file  is  created  on  FORTRAN  unit  10  and 
is  one  possible  input  to  VARVARYl.  FORTRAN  allows  usage  of  a  disk 
file  as  on  I/O  unit. 

(b)  A  file  of  the  corresponding  gradients  in  standard  for¬ 
mat  on  FORTRAN  unit  11. 

(c)  A  listing  containing  the  new  points  and  the  corre¬ 
sponding  program  values,  the  gradients,  and  the  computed  incre¬ 
ments  is  produced  on  the  printer. 

(d)  The  subsystem  includes  the  following  routines: 

1.  POINTCOMP. 

2.  PARTL. 

3.  PREPR. 

4.  The  program  to  be  tested. 

(2)  VARVARYl.  The  VARVARYl  component  of  the  POINTCOMP  sub¬ 
system  takes  the  file  created  by  POINTCOMP  on  unit  10  as  input  or 
a  similarly  structured  user  provided  file.  The  user  must  also  in¬ 
put  a  fila  containing  higher,  and  lower  bounds,  and  increments  .for.* 
each  variable.  VARVARYl  will  read  in  each  line  of  the  unit  10 
file  and  from  each  point  inputted  VARVARYl  creates  new  points  by 
varying  each  variable  in  turn  through  the  given  range  using  the 
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given  increment  keeping  the  other  values  fixed.  For  each  new 
point,  the  program  to  be  tested  is  evaluated  and  the  point  and  the 
corresponding  tested  program  value  are  printed  out. 

(3)  In  addition.  The  VARVARYl  routine  computes  the  gradient 
at  each  new  point  and  prints  its  components.  The  new  points  and 
the  corresponding  tested  routine  values  are  written  out  on  unit  12 
and  the  corresponding  gradients  are  written  out  on  unit  11.  When 
one  variable  is  varied  and  the  others  held  constant,  the  differ¬ 
ence  in  the  output  values  divided  by  the  difference  in  the  varied 
variable  values  is  called  the  difference  quotient.  As  the 
VAR\iARfl  routine  varies  a  variable  through  its  range,  various  com¬ 
binations  of  difference  quotients  are  computed  and  printed  for 
that  variable. 

(4)  The  output  from  the  VARVARYl  routine  includes: 

(a)  A  listing  of  the  new  input  points  and  the  associated 
output  values  of  the  routine  being  tested.  This  listing  also  con¬ 
tains  the  difference  quotients. 

(b)  A  listing  of  the  corresponding  gradients. 

(c)  A  file  containing  the  new  gradients  in  standard  fonnat 
on  unit  11. 

(d)  A  file  containing  the  new  generated  points  in  standard 
format  on  unit  12. 

(e)  A  scratch  file  containing  input  variable  values  and 
the  associated  output  values  on  unit  14. 

(f)  The  subsystem  includes  the  following  routines: 

1.  VARVARYl. 

2.  PARTL. 

3.  PRtPR. 

4.  The  program  being  tested. 

b.  ANALYZ  Subsystem.  This  subsystem  analyzes  a  file  of  gradi¬ 
ents  in  standard  format  on  unit  11  produced  by  the  POINTCOMP, 
VARVARYl,  or  GRID  subsystem.  The  statistics  produced  include: 

(1)  For  each  gradient,  the  ratio  of  every  component  to  the 
minimum  component  in  absolute  value. 


(2)  For  each  component  of  each  gradient  the  change  needed  in 
the  corresponding  variable  in  order  to  change  the  output  by  one 
unit. 

(3)  For  each  component  of  each  gradient: 

(a)  If  the  component  is  not  the  largest  in  absolute  value, 
the  ratio  of  the  component  to  the  absolute  value  of  the  maximal 
component. 

(b)  For  the  maximal  component,  the  ratio  of  the  absolute 
value  of  the  maximal  component  to  the  absolute  value  of  the  next 
largest  components. 

(4)  The  gradients  are  also  aggregated  component  by  component, 
and  statistics  (1)  to  (3)  described  above  are  computed  for  the 
aggregate  vector. 

(5)  The  output  from  ANALYZ  includes  a  listing  of  the  statis¬ 
tics  described  above. 

(6)  The  subsystem  includes  the  following  routine:  ANALYZ. 

c.  GRID  Subsystem.  This  subsystem  is  comprised  of  three 
linked  but  Independent  components— GRID,  REARRANGE,  and  DIFFQUOT. 

(1)  GRID.  The  GRID  subsystem  needs  a  range  and  step  size 
for  each  variable.  The  subsystem  varies  each  variable  from  the 
lower  bound  to  the  upper  bound,  incrementing  each  value  by  a  given 
step  size  to  obtain  the  next  value.,  thereby  producing  all  combinations 
of  variable  values.  The  routine  to  be  tested  is  evaluated  at  each 
point  and  the  gradient  is  computed  at  each  point  as  well. 

(a)  The  outputs  from  GRID  include: 

1_.  A  listing  of  the  new  points  and  the  associated  rou¬ 
tine  values. 

2^.  A  listing  of  the  gradients. 

2-  A  file  in  standard  format  containing  the  new  points 
and  the  associated  values  on  unit  10.  This  file  is  one  possible 
input  to  REARRANGE. 

4.  A  file  in  standard  format  containing  the  corresponding 
gradients  on  unit  11.  This  file  is  the  other  possible  input  to 
REARRANGE. 
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5.  A  scratch  file  containing  the  new  points,  the  asso¬ 
ciated  vaTues  and  some  pointers  on  unit  12.  This  file  is  used  as 
an  input  to  REARRANGE. 

6.  A  scratch  file  containing  the  gradients  and  tested 
routine  values  on  unit  13. 

(b)  The  components  of  this  subsystem  are: 

1.  GRID. 

2.  PARTL. 

3.  PREPR. 

4.  The  program  to  be  tested. 

(2)  REARRANGE.  This  routine  reads  from  unit  14,  one  of  two 
files  produced  by  GRID  on  units  12  or  13.  The  routine  also  util¬ 
izes  the  file  produced  by  GRID  on  unit  12.  The  output  data  is  re¬ 
arranged  variable  by  variable  to  facilitate  further  analyses. 

(a)  The  output  from  REARRANGE  includes: 

A  listing  of  the  rearranged  input  file,  where  adjoin¬ 
ing  lines  in  the  same  section  differ  only  in  one  variable  value, 
each  point  is  numbered  to  facilitate  linkage  of  this  output  to  the 
GRID  output. 

A  file  containing  similar  information  produced  on 

unit  15. 

(b)  The  subsystem  has  one  component:  REARRANGE. 

(3)  DIFFQUOT.  This  routine  reads  the  GRID  output  files  pro¬ 
duced  on  units  10  and  12.  The  routine  outputs  a  list  of  differ¬ 
ence  quotients  for  the  variable  requested.  The  outputted 
quotients  are  headed  by  the  numbers  of  the  relevant  GRID  points, 
facilitating  the  linkage  between  the  GRID  output  and  this  output. 
The  subsystem  has  one  component:  DIFFQUOT. 
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1-2.  DESCRIPTION  OF  SUBSYSTEMS.  Subsystem  components  are  shown 
in  Table  1-1. 


Table  1-1.  Subsystem  Composition 


Subsystem 


Major 

components 


Minor 

components 


POINTCOMP 

POINTCOMP 

PARTL 

VARVARYl 

PREPR 

ANALY? 

ANALY Z 

— 

GRID 

GRID 

PARTL 

REARRANGE 

DIFFQUOT 

PREPR 

1-3.  SUBSYSTEM  INPUTS  AND  OUTPUTS.  Figure  1-1  and  Table  1-2  show 
inputs,  outputs,  and  intercommunication  between  the  subsystems. 


PARTYl.  — - 
PRtPR 


POINTCOMP 


GRID 


■  PARIYl 
.  PREPR 


:  Units  10,  12 


Unit  11 

( 

1 

1 

■  Unit  U 

o' 

j  Units  10,  12 

/ 

1 

, 

PAKfYL  •  -  »  /A 

V.ARVARYl 

ANAEYZ 

1'  i 

REARRANGE 

I 

f 

DIFFQUOT  ■= - *■  PREPR 

PREPR 


UOTES:  POINTCOMP.  VARVARYl,  and  GRID  routines  each  may  be  run  independently  of  the  other 
major  routines.  The  other  Xalor  routines  may  be  run  independently  once  POINTCOMP  or  GRID  have 
been  run  tn  get  up  the  required  input  tiles,  as  indicated  above. 

The  head  of  the  arruH  denotes  input;  the  tail  of  the  arrow,  output. 


Figure  1-1.  Data  Flow 


1-5 


Table  1-2.  Unit  Utilization  (all  routines  utilize  Units  5  and  6 
in  addition  to  those  listed  below.) 


Routine 

name 

Unit 

number 

Type 

Format 

Description 

POINTCOMP 

10 

Output 

Standard 

Each  line  contains  input  variable  values 
and  the  corresponding  value. 

11 

Output 

Standard 

Each  line  contains  the  corresponding  gradients. 

VARVARYl 

10 

Input 

Standard 

File  created  by  PONTCOMP,  containing  the 
the  variable  values. 

11 

Output 

Standard 

Gradients  corresponding  to  new  points 
generated  by  VARVARYl. 

12 

Output 

Standard 

New  points  generated  by  VARVARYl  and  the 
associated  value  of  the  tested  routine. 

13 

Output 

Nonstandard 

Scratch  file  containing  gradients. 

14 

Output 

Nonstandard 

Scratch  file  containing  input  variable  values. 

ANALYZ 

11 

Input 

Standard 

Gradient  file  created  by  POINTCOMP,  VARVARYl, 
or  GRID 

GRID 

10 

Output 

Standard 

File  containing  input  variable  values  and 
the  associated  values  of  the  tested  routine. 

11 

Output 

Standard 

File  containing  the  corresponding  gradients. 

12 

Output 

Nonstandard 

Scratch  file  containing  the  input  variable 
values  tested  routine  value  and  a  pointer  to 
the  variable  varied  to  create  the  point. 

Used  as  input  to  REARRANGE. 

13 

Output 

Nonstandard 

Scratch  file  containing  the  corresponding 
gradients  and  the  value  of  the  tested 
routine.  Used  as  input  to  REARRANGE. 

REARRANGE 

10 

Input 

Standard 

File  produced  by  GRID,  containing  the 
input  variable  values  and  tested  routine 
value  on  Unit  12,  containing  the  indication 
of  which  variables  were  varied. 

12 

Input 

Nonstandard 

File  produced  by  GRID  on  Unit  12,  containing 
the  indication  of  which  variable  was  varied. 

Used  in  the  reordering  computations  by  REARRANGE. 

14 

Input 

Nonstandard 

The  file  to  be  reordered,  this  file  could 
be  the  file  produced  by  GRID  on  Units 

12  or  13. 

15 

Output 

Nonstandard 

A  file  containing  the  rearranged  output. 

OIFFQUOT 

10 

Input 

Standard 

File  containing  input  variable  values, 
produced’ by  GRID.* 

12 

Input 

Nonstandard 

File  produced  by  GRID,  indicating  which  variables 

were  varied  to  create  each  point. 
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The  following  (Figure  1-2)  outlines  the  units  utilized  by  the  dif¬ 
ferent  routines. 


NOTE:  Numbers  denote  I/O  unit  numbers;  primes  denote  same  unit 
but  different  formats. 


Figure  1-2.  Subsystem  Unit  Utilization 


1-4.  STANDARD  FILE  STRUCTURE.  A  standard  file  has  the  following 
structure: 

a.  Line  1  is  Comprised  of  Fields  Al,  A2,  and  A3 

(1)  Al  is  the  number  of  points  or  gradients  in  15  format. 

(2)  A2  is  the  number  of  variables  in  13  format. 

(3)  A3  is  a  positive  number  in  F10.7  format.  Any  vector 

components  smaller  than  this  number  in  absolute  value  will  be  con¬ 
sidered  to  be  zero. 

b.  Line  2  (format).  This  line  indicates  the  format  that  was 
used  to  write  each  data  line  in  the  file  and  which  may  be  used  to 
read  the  subsequent  lines  in  the  file. 

c.  Other  Lines.  These  lines  contain  data  consisting  either  of 
variable  value  vectors  (points)  or  gradients.  The  leftmost  column 
gives  data  on  variable  1,  the  next  column  to  the  right  on  variable 
2,  etc. 
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CHAPTER  2 


THE  POINTCOMP  SUBSYSTEM 


Section  I.  THE  POINTCOMP  ROUTINE 


2-1.  INTKUUUCTIUN.  This  routine  uses  the  gradient  to  pick  new 
points  which  miniinize/inaxiinize  (locally)  the  value  of  the  output 
of  the  routine  being  tested.  In  addition,  the  gradient  at  each 
point  is  also  printed  out.  The  analyst  may  use  the  output  to 
check  out  tiie  routine  being  tested  in  the  following  manner: 

a.  The  analyst  can  examine  each  variable  to  see  if  it  must  be 
i ncreased/ decreased  to  increase/decrease  the  output.  The  analyst 
can  also  check  for  attainment  of  values  from  which  no  increase/ 
decrease  of  the  output  variable  is  possible.  These  results  should 
be  explainable  intuitively. 

b.  The  analyst  may  examine  each  gradient  to  determine  the 
relative  effect  of  each  variable  at  that  point. 

c.  The  analyst  may  examine  the  components  of  the  gradients  for 
all  points  to  analyze  the  marginal  returns. 

d.  Changes  in  the  sign  of  a  component  of  a  gradient  from  one 
point  to  the  next  indicate  that  a  potential  local  optimum  for  that 
variable  lies  between  the  two  successive  values  of  the  coordinate. 

e.  POINTCOMP  calls  PARTL  to  compute  partial  derivatives  nu¬ 
merically.  POINTCOMP  calls  PREPR,  a  user  provided  driver  routine 
which  obtains  an  output  value  from  the  program  being  tested. 

2-2.  BACKGROUND 


If  •••.  x^)  is  a  real  valued  function  of  several  variables, 
•••.  Xj,)  =  (■g^^(xi,  ...»  x„),  ...,|l(xi,  ...,  Xp)) 

points  in  the  direction  of  greatest  local  increase  in  f  from  (x,, 
...,  x^),  and  -yf(xp  ...,  x  )  points  in  the  direction  of  great¬ 
est  local  decrease  in  f  from  (Xj^,  ...,  Xp).  We  justify  this  fact 
by  the  following  argument: 


If  (y^ . Yp)  is  a  point  "close"  to  (xj . Xp), 

f(yi . yp)  -  f(xi . Xp)  is  approximated  by: 

H  (xi.  ....  x„)  (Yi  -  x^) 
i=l  ^^i 

=Vf(x^ . Xp)  •  (y  -  x) 

=  IIVf(xi,  ....  Xp)i|  ||y  -  x||cos  e 
where:  y  =  (yi*  ....  yp); 

X  “  (Xj.  ...»  Xp); 

|(A(|  is  the  magnitude  of  a  vector  A;  and 
0  is  the  angle  betweenVf(xj,  ...»  Xp)  and  (y  -  x) . 
So 

f(yj^.  ...»  Yp)  "  f(xj^»  ...»  Xp)  ||Vf(xj^.  ...»  Xp)||  I  jy 
and  therefore: 


-  xl I  cos  0, 


f(yii 


yp)  -  f(xi»  ...»  X  ) 

Tj - r -  =llVf(xi,  ..•»  X  )||cos  0 

l|y  -  x||  ^ 


Now  -1  i  cos  0  i  1,  so 

max  f(yi  •••»  yri)  ~  f^^i  ...»  Xp) 

I ly-xlij'u  My  -  x|| 


1 1  Vf  •  •  •  >  ’^nM  I  • 


Therefore  the  greatest  change  in  the  function  f  per  unit  distance 
between  y  and  x  (i.e.»  | |y  -  x||)  is  the  magnitude  of  the  gradient 
iivn  Xi ,  . . . ,  Xp) 1 1 .  This  maximum  is  reached  for  cos  0  =  1  or  0 
=  0,  soothe  direction  of  greatest  increase  is  0  =  0»  so  y  -  x  and 

f(Xi»  ...»  Xp)  are  col  linear.  Therefore,  Vf(xi . Xp)  gives 

both  the  magnitude  and  direction  of  the  maximal  change. 


L i kewise: 


min  f(yi  •••»yn)  "  f^^l  ^^n) 

My-xM?'o  My  -  ’^M 


"llVf(X]^  ...»  Xpll 
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In  this  case,  cos  0  =  -1,  so  0  =  180“  and  the  direction  of  great¬ 
est  decrease  (least  increase)  per  unit  distance  between  y  and  x 
dly  -  xl|)  is  -Vf(xp  Xj^),  the  magnitude  of  greatest  de¬ 
crease  is  llVf(xi.  .i.,  x„)i|. 

2-3.  DISCUSSION  OF  METHODOLOGY.  Several  criteria  are  used  in  the 
program  to  decide  upon  the  distance  to  be  moved  along  the  gradient 
in  order  to  obtain  the  next  point.  These  criteria  are: 

a.  The  maximum  change  in  one  step  for  each  variable  (an  input 
parameter)  divided  by  the  size  of  the  component  of  the  gradient 
corresponding  to  that  variable. 

b.  First,  compute  a  certain  fraction  of  the  difference  between 
the  current  variable  value  and  the  bound  on  the  variable  in  the 
appropriate  direction  (an  input  parameter)  whenever  the  bound  ex¬ 
ists.  When  the  variable  isn't  bounded,  use  twice  the  current 
value  as  a  bound.  Next,  the  number  just  described  is  divided  by 
the  size  of  the  corresponding  component  of  the  gradient. 

c.  Compute  the  minimum  of  the  numbers  developed  in  parts  a  and 
b  and  the  default  multiplier  read  in  as  an  input  parameter  on  the 

8th  input  card  (see  pg  2-6).  This  step  computes  a  gradient  multiplier. 

d.  A  candidate  increment  is  determined  by  multiplying  the  gra¬ 
dient  by  the  number  derived  by  the  process  described  in  c.  If  the 
gradient  multiplier  is  too  small,  use  the  A7  field  on  card  1.  If 
necessary,  the  increment  is  multiplied  by  a  number  sufficient  to 
increase  each  component  to  the  minimum  threshold.  These  minimum 
thresholds  are  also  input. 

e.  Lastly,  the  increment  vector  is  added  to  the  original 
point,  obtaining  a  new  point  candidate.  The  components  of  this 
new  point  exceeding  the  bounds  of  the  corresponding  variable  (if 
any)  are  set  equal  to  the  value  of  the  bound,  and  the  result  is 
taken  to  be  the  next  point. 

2-4.  LIMITATIONS.  As  currently  compiled,  the  POINTCOMP  routine 
is  subject  to  the  following  limitations: 

a.  Testing  of  the  program  must  proceed  by  testing  one  out¬ 
put  variable  at  a  time  if  the  program  has  several  output  vari¬ 
ables. 

b.  If  any  variables  input  to  the  program  being  tested  are 
integers,  it  is  possible  that  their  incrementation  and  effect  upon 
the  output  variables  will  be  reduced  due  to  truncation. 
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c.  Current  limits  are  set  to  20  variables  and  500  generated 
points  at  most.  This  is  easily  changed. 

d.  Testing  Monte  Carlo  programs  is  feasible,  but  potentially 
time  consuming.  The  PREPR  routine,  called  by  POINTCOMP,  must  call 
the  routine  being  tested  repetitively  and  average  the  output  val¬ 
ues  before  returning.  Replication  is  necessary  in  order  to  aver¬ 
age  out  random  effects.  A  better  approach  is  to  use  the  Monte 
Carlo  variables  as  ordinary  input  variables  and  let  POINTCOMP  as¬ 
sign  their  values  with  no  replication. 

e.  The  function  represented  by  the  routine  to  be  tested  should 
be  well-behaved  (i.e.,  differentiable). 

f.  The  order  in  which  the  criteria  for  step  size  are  examined 
is  fixed,  so  there  is  a  hard  wired  priority  among  the  criteria. 

2-5.  RUN  SETUPS 

a.  Developing  an  Absolute  File  in  ASCII 


@MAP,S  .Name  of  absolute  element 

IN  03PR0GTEST.  POINTCOMP 

IN  O3PR0GTEST.  PARTIAL 

IN  Element  containing  PREPR  subroutine 

IN  Programs  to  be  tested 

LIB  LIB$*FTN3. 

END 

b.  To  Execute 


0USE  10,  file  name  to  contain  variable  values  and 

values  of  the  tested  routine 

@USE  11,  data  file  name  for  gradient  file  in  standard  format 

@ASG,A  First  file  name 

@ASG,A  Second  file  name 

@XQT  Absolute  element  created  by  the  PMAP 

Input 

Deck 


c.  Input  Deck.  The  following  tables  describe  the  input  deck 
(Table  2-1),  show  a  sample  input  deck  (Table  2-2),  and  a  sample 
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Table  2-1.  Description  of  Input  Deck 
(page  1  of  2  pages) 


Line  1 


Field:  A1  A2  A3  A4  A5  A6  A7  A8 


Format:  13  13  12  F10.7  F10.7  F12.0  F12.5  13 


where  the  fields  are  defined  as  follows: 

A1  is  the  maximum  number  of  points  to  be  generated. 

A2  is  the  number  of  variables. 

A3  is  +1  if  output  is  to  be  increased,  -1  if  output  is  to  be 
decreased. 

A4  is  a  positive  real  number.  Any  number  smaller  than  this 
number  will  be  considered  to  be  zero. 

A5  is  a  positive  real  number,  gradients  smaller  than 
this  number  will  be  considered  to  be  zero. 

A6  is  a  code  number.  This  number  can  be  used  to  indicate 
that  a  variable  is  unbounded  above  or  below. 

A7  describes  the  minimal  multiplier  desired  for  the  gradi¬ 
ent. 

A8  the  run  is  executed  in  a  debugging  mode  (more  printout) 
when  this  field  contains  a  -1. 


Line  2.  (Format  for  reading  in  each  of  the  following  six  lines, 
lines  3-8  which  follow.  Parentheses  must  bound  the  format 
as  shown.) 


Line  3.  Initial  values,  one  for  each  variable,  in  the  format  de¬ 
scribed  by  line  2. 
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Table  2-3.  Description  of  Input  Deck 
(page  2  of  2  pages) 


Line  4.  Lower  bounds  allowed  for  each  variable,  in  the  format  de¬ 
scribed  by  line  2.  If  some  variable  is  unbounded  below, 
indicate  this  by  using  the  number  from  field  A6  of  the 
first  line. 


Line  5.  Upper  bounds  allowed  for  each  variable.  Rest  of  descrip¬ 
tion  as  in  line  4. 


Line  6.  The  preferred  maximum  single  step  change  in  each  vari¬ 
able,  one  number  for  each  variable. 


Line  7.  The  preferred  minimum  single  step  change  in  each  vari¬ 
able,  one  for  each  variable. 


Line  8.  Default  gradient  multiplier  for  each  variable,  one  number 
for  each  variable.  These  will  only  be  used  for  unbounded 
variables  and  are  only  useful  to  reduce  the  number  by  which 
the  gradient  will  be  multiplied. 


Line  9.  (Format  for  writing  one  line  of  variable  values  and  the 
corresponding  program  output  value  on  unit  10.  This  fonnat 
should  accommodate  at  least  n+3  values,  where  n  is  the  num¬ 
ber  of  variables.) 

Line  10.  (Format  for  writing  out  one  line  of  variable  values  and 
the  corresponding  program  output  value  on  the  printer.) 
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Table  2-2.  Sample  Input  Data 

4^  9^+1^  ,001^^  .001®  1001.^  .1000^+1^ 

(13F5.0)  format  for  inputting  the  following  six  lines 


5.0 

0.5 

0.5 

0.5 

Initial  value 

0.0 

0.0 

0.0 

0.0 

Lower  bounds 

1000.0 

1000.0 

1000.0 

0.9 

Upper  bounds 

10.0 

10.0 

10.0 

0.5 

Preferred  maximum  step 

0.0 

n.o 

0.0 

n.o 

Preferred  minimum  step 

1.0 

1.0 

1.0 

1.0 

Default  multiplier 

(not  really  needed  for 
bounded  variables) 


First  column  describes  variable  1. 

Second  column  describes  variable  2. 

Etc. 

(6(7X,F13.5)) 

(lX,6(fX,F13.5)) 

^Number  of  variables. 

^Maximum  number  of  points. 

^Increase  the  output. 

^Numbers  lower  than  this  are  considered  zero. 
^Gradients  smaller  than  this  are  considered  zero. 
^Code  number  indicating  no  bounds. 

^Minimal  gradient  multiplier. 

^Not  in  debug  mode. 


Format  for  writing  values  and 
output  on  unit  10 

Format  for  writing  values 
and  output  on  the  printer 
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Table  2-3.  Sample  Run  Deck 


PUSH  10,03MAT1. 

PASG.A  03MAT1. 

PUSE  11,03MAT2. 

PASG.A  03MAT2. 

PXQT  absolute  deck  created  by  PMAP 


mm  +1 

(13F6.0) 

.001 

5.0 

0.5 

0.0 

0.0 

1000.0 

1000.0 

10.0 

10.0 

0.0 

0.0 

1.0 

1.0 

(6(7X,F13.5)) 

(1X,6(7X,F13.5)) 


.001 

1001. 

0.5 

0.5 

0.0 

0.0 

1000.0 

0.9 

10.0 

0.5 

0.0 

0.0 

1.0 

1.01 

.1000  +1 


2-6.  OUTPUT  DESCRIPTIONS  AND  SAMPLE  OUTPUT 

a.  Tuples.  Each  row  corresponds  to  a  point  where  the  right¬ 
most  column  is  the  value  of  the  routine  being  tested  and  the  left¬ 
most  column  represents  variable  1,  the  column  to  its  right  vari¬ 
able  2,  etc. 

b.  Gradients.  Each  row  is  the  gradient  at  the  point  described 
by  the  corresponding  row  of  the  tuples  output,  e.g.,  the  first  row 
is  the  gradient  at  the  point  described  by  the  first  row  of  the 
tuples  printout. 

c.  Increments.  Each  row  is  the  computed  vector  difference  of 
the  corresponding  tuples  (except  where  boundaries  are 
encountered).  The  first  row  is  the  second  tuple  minus  the  first, 
the  second  row  is  the  third  tuple  minus  the  second,  etc.  The  fol¬ 
lowing  (Table  2-4)  exhibits  sample  printer  output. 
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Table  2-4.  Sample  Output  to  the  Printer 


Tuples 


5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

.61902 

.75809 

.63333 

.99916 

5.00119 

.73820 

.97516 

.76667 

1.24159 

5.00270 

.80924 

1.08781 

.84667 

1.37902 

5.00270 

.84061 

1.13401 

.88222 

1.43833 

5.00270 

.85178 

1.14995 

.89492 

1.45926 

5.00270 

.85513 

1.15468 

.89873 

1.46551 

5.00270 

.85513 

1.15590 

.89873 

1.46626 

5.00270 

.85513 

1.15715 

.89975 

1.46753 

Gradi ents 

.00048 

.30919 

.67048 

.34638 

.00359 

.35940 

.65458 

.40207 

.00852 

.40086 

.63565 

.45144 

.01236 

.42295 

.62294 

.47942 

.01428 

.43221 

.61698 

.49139 

.01500 

.43545 

.61481 

.49560 

.01522 

.43641 

.61415 

.49686 

.01523 

.43631 

.61415 

.49674 

.01527 

.43670 

.61405 

.49663 

Increments 

.00000 

.11902 

.25809 

.13333 

.00119 

.11919 

.21707 

.13333 

.00151 

.07104 

.11264 

.08000 

.00000 

.03137 

.04620 

.03556 

.00000 

.01117 

.01594 

.01270 

.00000 

.00335 

.00473 

.00381 

.00000 

.00000 

.00122 

.00000 

.00000 

.00000 

.00126 

.00102 
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d.  Output  File  Descriptions 
(1)  File  on  Unit  10 


(a)  Line  1 

Fields 

A1 

A2 

A3 

Formats 

15 

13 

F10.7 

where 

A1  is  the  number  of  points  in  the  file. 

A2  is  the  number  of  variables. 

A3  is  the  number  determining  when  small  numbers  are  considered 

to  be  zero. 

(b)  Line  2.  (Format  for  reading  in  one  point  and  its  cor¬ 
responding  program  value.) 

(c)  Other  Lines.  Each  line  contains  a  set  of  variable 
values  and  tbe  corresponding  tested  program  value. 

(2)  File  on  Unit  11 

(a)  Same  as  file  on  unit  10. 

(b)  Same  as  file  on  unit  10. 

(c)  Other  lines  -  Each  line  contains  the  corresponding 
gradient. 
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e.  Sample  Output  File — Unit  10 


9  4  .0010000 


(6(7X,F13.5)) 

Tuples 

5.00000 

.50000 

.50000 

.50000 

5.00000 

.61902 

.75809 

.63333 

5.00119 

.73820 

.97516 

.76667 

5.00270 

.80924 

1.08781 

.84667 

5.00270 

.84061 

1.13401 

.88222 

5.00270 

.85178 

1.14995 

.89492 

5.00270 

.85513 

1.15468 

.89873 

5.00270 

.85513 

1.15590 

.89873 

5.00270 

.85513 

1.15715 

.89975 

f.  Sample 

Output  File- 

-Unit  11 

9  4  .0010000 

(6(7X,F13.5)) 

.00048 

.30919 

.67048 

.34638 

.00359 

.35940 

.65458 

.40207 

.00852 

.40086 

.63565 

.45144 

.01236 

.42295 

.62294 

.47942 

.01428 

.43221 

.61698 

.49139 

.01500 

.43545 

.61481 

.49560 

.01522 

.43641 

.61415 

.49686 

.01523 

.43631 

.61415 

.49674 

.01527 

.43670 

.61405 

.49663 

.73821 

.99916 

1.24159 

1.37902 

1.43833 

1.45926 

1.46551 

1.46626 

1.46753 
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2-7.  POINTCOMP  ROUTINE  LISTING 


PARAMETER  F ILEl^lQ ,F ILE2:ll.N0VALS:20tN0PTS  =  500 i INFILE=5 1 OUTFILrS 
PARAMETER  NOVALI=NO VAtS«l 
C 

CHARACTER *80  FORM Tl# F0RMT2tF0RMI3 

C 

DIMENSION  FSIVALtNOVALSIfLOUBOStNOVALSI tHYBDS INOVALS  I 
DIMENSION  MAXCHGINOVALSI tCRADINOPTStNOv ALS>t TUPLES INOPTSfNOVALl I 
DIMENSION  0ELJAl(N0VALS)t0ELTA2<N0VALSI iDELTAINOPTStNCVALSI 
DIMENSION  SIGNINOVALSI tINPUT INOVALS  It  MIN  INOVALS )» NOMULT INOVAL SI 
C 

INTEGER  NOVAHStUPOOUNtMAXPTStSIGNtltJtPfB. SWITCH 
C 

REAL  EPSLONtFLATtNOLTM tFSTVALtLOWBDStHVBDS  tMAXCHB 
REAL  GRAOtSUMtiNPUTt INCRtPARTTL t DEL T Alt  DEL TA2f DELTA 
REAL  TUPLEStVALUE tNOMULTtMINtMINltNINMLTtMAX 
C 

REWIND  FILEl 
REWIND  FILE2 

REA 01  infile tLOOOO I  NOV ARS t MAXPTS t UPDOWN t EPSLON tFL AT  tNOLIMt 
1  MiNMLTtSWirCH 

10000  F0RMATI2I3iI2t2Fia.7  t2F12.0tI3l 


IF  {SWITCH  .NE.  -II  GOTO  30010 

WRITEI0UTFILt200S0l  NO VARS t MAXP TS tUPDOWN tEPSLON tFL AT t NOL IM,MINML T 
20060  FORMAT!’  *»IM V ALS • t/ t IX 1 213 1 12 t2FlD. 7 1 F12.0 t F12 . 5 1 
30010  CONTINUE 

C  •**•* 

READIINFiLEiinOlOl  FORMTl 
10010  FORMATI A80I 

READ! INFiLEt FORMTl I  IFSTVALI 1 1 tl-ltNOVARS» 

REAOIINFIlE  t FORMTl I  IL OWBDSI 1 1 1 1=1 1 NO VARSI 
READ! infile t FORMTl I  IH YBOSI I I t 1= It NOV ARS » 

REAOIINFILE  t FORMTl I  I H AXCHG 1 1 1 1 1  =  1 1 NO V ARS) 

READI INFILF  t FORMTl I  I M IN  1 1 1 1 1= 1 1 NO V ARS I 
REAOIINFILE .FORMTl  I  INOMULTI 1 1 1 1  =  1 . NO VARS I 
REAOIINFILE. lODlOJ  F0RMT2 
READI INFIlE .100101  F0RMT3 
C 

n  ISWITCH  .NE.  -II  GOTO  3002Q 
WRI rEIOUTFIL.210001  IFSTVALI 1 1 .1=1. NO VARSI 
21000  FORMATI*  * • FS T V AL * « ’ tlUl / . IX .6F20 . 10  1 1 

MRITEI0UTFIL.21010I  ILOWBDSI I ) .1=1. NO VARSI 
21010  FORMATI’  •.LaVBDS • • *  .101 / . IX .6F20. 10 1 1 

WRI TE I0UTFIL.21O20I  IHYBOSIII .1= l.NOVARS  I 
21020  FORMATI*  ••HY«0S***.10I/.1X.6F20.10I  I 

WRITCIOUTFrL.210301  IN AXCHGI I I . 1=1 . NC VARSI 
1030  FORMATI’  *  •  MITXCHG  •  •  ’  .1 0 1 /.  IX  .6F  20  .  ID  1 1 
WRI TE I0UTFIL.2U06SI  I M IN  1 1 1 • 1= 1. NOV ARS I 
20065  FORMATI*  • * MTNS* « ’ . 1 0 I /. 1 X .6F 2 0. 10 I i 

WRI TE I0UTFIL.2U067  I  INOMULTI I ) .1  =  1 . NO VARS I 
20067  FORMATI’  •• NOMULT • « «  .101 / . IX .6F20.1D I  I 
30020  CONTINUr 
c  •••*• 
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C  INITIALIZATION 

DO  100  I=1>N0VARS 
TUPLESIltII=r.STVALm 
INPUT«II=rST»AHI  » 

100  CONTINUE 

CALL  PREPR(INPUTtVALUE) 

TUPLES* It  NOV ALII zVALUc 
C  LOOP  TO  CRfATE  POINTS 

DO  3000  P=2flfAXPTS 
DO  300  I=ltNOVARS 
INPUTIII-TUPLESIP-lt II 
300  CONTINUE 

DO  ADO  I=ltNCVARS 

CALL  PARTL*  NOVARStl. INPUTtEPSLONf PARTTLt INCRiTUPLESfP-ltNOVALlll 
C 

IF  (SWITCH  .HE.  -II  GOTO  30030 

WRITE (OUTFIL *200001  I INPUT ( J 1 1 J= ItNOV ARS 1 1 P ART YL i INCR  t 
1  TUPLESIP-ltNOVALl  I 

20000  FORHATI//.*  ••INPUT  PARTIAL^^ • »10 (/ tlX . GF2D. ID  1 1 
30030  CONTINUE 

C  ••••• 

GRADIP-ltII=P«RTyL 
AOO  CONTINUE 

00  500  I=ltNftVARS 
SI6NIII=UPOOUN 

IF  IGRACIP-ltll  .LT.  0.1  SIGN( II =-UPDOWN 
500  CONTINUE 

C  ••••• 

IF  (SWITCH  .NE.  -II  GOTO  300AO 
HRITE(0UTFIL.20010l  (GRAD ( P-l.II . 1=1 t NO V ARS I 
20010  FORMAT*//. •  ••GRA0*» • . 10 ( / t IX »6F20.10) I 
HRITE(0UTFILt2D015l  ( SIGN ( 1 1 • 1=1 .NO VARS  I 
20015  FORMAT*//. •  ••SIGN** ’. 10 ( /.IX . 101311 
300AO  CONTINUE 
C  ••••• 

C  COMPUTE  CANDIDATE  MULTIPLIER  BASED  MAX  CHANGE  PER  COORDINATE 
DO  600  I=1.N0VARS 
OELTAK  I1  =  ABS'IN0MULT*II> 

IF  (AeS(GRA0(P-l,I|l  .LT.  EPSLONI  GOTO  600 
OELTAK It  =  AeS(MAXCH6 (I  1 1 /ABS ( GRA0( P-1 . 1 1  I 
600  CONTINUE 

C  ••••• 

IF  (SWITCH  .NE.  -II  GOTO  30050 
HRITE(0UTFIL.20a2UI  (OELT A1 ( I I ,1=1 .NO VARSI 
20020  FCRHATf//.*  *»DEL TAl • • *.10 ( / .iX.6F20.10 I  I 
30050  CONTINUE 
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c  *«••• 

C  COMPUTE  CANI>IOATE  HULTIPLIER  BASED  ON  DISTANCE  FROM  BOUNDS 
C  OR  DOUBLING  THE  VALUE  IF  NO  BOUND  EXISTS 

DO  800  I=liNOVARS 

DELTAZtllrAMlNKABSt  TUPLESIP-1«T)/ABSIDRADIP-1»IM  IfNOHULTlIM 
IF  (SI6NIII  4.T.  0  I  GOTO  700 

IF  (ABSf HYBOStll-NOLIMl  .GE.  EPSLON)  DELTA2(I)= 

1  HFLOATIPI-I. 1/IFL0AT«P|41. 1 1*1 ABSIHTBOSi l)-TUPLESIP-l*It 11 

2  /ABSIORAO(P-ltl)l 
GOTO  800 

700  IF  (ABSILOUBDSdI-NOLIMI  .GE.  EPSLONi  DELTAZdlr 

1  KFLOAUPI-I;.  I/<FL0AT(PI>1.  n*t  ABS(LCU8DSd»-TUPLESIP-ltI»>  ) 

2  /ABSiGftADIP-l.Di 
800  CONTINUE 

C  ••••• 

IF  (SWITCH  .NE.  -1>  GOTO  30060 
URirE(0UTFIL<2004GI  (0ELTA2d  )  d=l  >NCVARSt 
200*6  FOfiMAT(//t‘  <»DELTA2**».10l/»lXt*F30.10> > 

30060  CONTINUE 

C  •••*• 

C  COMPUTE  THE  INCREMENTATION  VECTOR 

00  900  IrltNCVARS 

DELTA(P-l.Il=AMINl(OELTAldltOELTA2(H  J 
900  CONTINUE 

C  •••«* 

IF  (SWITCH  .NE.  -i»  GOTO  30070 
HRITE(0UTFlLte0045l  (DELT A(P-l »I ) »I=1 tNO VARS » 

200*5  FORMATf//**  •ADELT A. 0* • 10( / 1 1X»6F20.10 ) I 
30070  continue 

C  •••«• 

C  TO  COMPOTE  THE  MINIMUM  GRADIENT  MULTIPLIER 
MIN1=0ELTA(P-1,1) 

00  925  I:l»NO!VARS 

IF  (0ELTA(P-1jH  .LT.  EPSLONI  GOTO  925 
IF  (0ELTA(P-Idl  .LT.  MINII  MINI  =  0ELTAIP-1  »I  I 
9Z5  CONTINUE 
C 

IF  (SWITCH  .NE.  -11  GOTO  30080 

WRITE(0UTEIL. 200**1  (DELTA (P-I . I ) » 1  =  1  .NO VARS  I  .MINI 
ZOO**  FORMATC/'/.*  ••DELTA. 1***.10(/.1X.6F20. 101) 

30080  CONTINUE 

C  TO  COMPUTE  A  CANDIDATE  INCREMENT 

IF  (MINI  .LT.  EPSLONI  MINlrMINMLT 
DO  950  I=1.N0VARS 

DELT*(P-1.II=HIN1^ABSIGRA0IP-1.I|  I 
950  CONTINUE 

C  ••••• 


>'-I4 


IF  (SWITCH  .ME.  -11  GOTO  30090 

HRITftOUTFILt 200401  lOELT* (P-1  til 1 1=1 tNO VARS ) »NIN1 
20040  FORHAT(//«*  ••CELT A. 2 1 HINl** • • 10 1 /t IX tEFZO .lU t  I 
30090  CONTINUE 

C  ••••• 

C  TO  GUARANTEE  A  HINlMAL  STEP  IN  EACH  COORDINATE 
HAX=0. 

DO  990  I=liN»VARS 
C  SET  MINISCULE  STEPS  TO  ZERO 

IF  IDELTA(P-ltl)  .GE.  EPSLONI  GOTO  9S0 
DELTA(P-ltI|:a. 

GOTO  990 

C  CHECK  TO  SEE  IF  CANDIDATE  INCREMENT  EXCEEDS  MINIMUM 
380  IF  (  OELTAIP-liIl  .GE.  MINIII  I  GOTO  990 

C  ••••• 

IF  (SWITCH  .KE.  -It  GOTO  30200 

WRITE  (0UTFILf‘2020ai  It  MAX  tMIN(  1 1  /  ABSI  DELTA  (P-1 1 1 1 1 
202U0  FORMAT(/flOO(/t*  •*NAX ••  •  tIS 1 4F20.10  )  ) 

30200  CONTINUE 

C  •••*• 

IF  ((MlN(Ii/OELTA(P-l.II>  .LE.  MAX)  GOTO  990 
MAX=MIN(I)/0EX.TA(P-1  .1 1 
990  CONTINUE 

IF  (  max  .LT.  EPSLONt  GOTO  975 
C  INCREASE  THE  STEP  SIZE 
DO  970  I=liNeVARS 
0ElTA(P-1.H=MAX«0£LTA(P-1»I) 

970  CONTINUE 
C  ••••• 

IF  (SWITCH  .nE.  -II  GOTO  30110 

NRITE(0UTFILt20042l  (CELT A (P-1 tl } t 1=1 tNOVA RS I t M AX 

20042  FORNATI/Zt’  *»DELT A. 3« M AX 1 10( / t IX t GF20. 10 1 1 
30110  CONTINUE 

C  •••*• 

37S  DO  985  I=ltN»VARS 

IF  (  DELTA(P-l.I)  .GT,  MAXCHGdI  I  DELTA  (P-1 1 1  l=MAXCHO  (I  I 
0ELTA(P-lfII=0ELTA(P-l,I)*SI6N(I» 

385  CONTINUE 

C  ••••• 

IF  (SWITCH  .ME.  -It  GOTO  30120 
WRITE(OUTFILrZ0043l  IOELTA(P-ltII .I=ltNOVARSI 

20043  FORNAT(//t’  ••DELTA.4*** tl0(/tlXt$F20.10n 

30120  CONTINUE 

C  ••••• 

C  CHECK  tuple  VALUES  AGAINST  THE  BOUNDS 
C  IF  ANY 

DO  1100  I=ltlU)VARS 
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TUPLESIPtH  =  TUPLES<P-l.II*OELT*«P-l.II 
IF  (ABSIHYBOSIII-NOLIN)  .LT.  EPSLONI  GOTO  1000 
IF  ITUPLESIPtI)  .GT.  HY8DSII))  TUPLES ( P t II =HYBDS 1 1 > 

1000  IF  lABSlLCUBOiSdl-NOLIKt  .LT.  EPSLONI  GOTO  1100 

IF  (TUPLES(P>JI  .LT.  LOUBOSdll  TUPLES  I P  >I  l=LOUBOS  II I 
1100  CONTINUE 
C  «•••• 

IF  ISUITCH  .NE.  -II  GOTO  30140 
URITEI0UTF1L*2004TI  ITUPLCSIP t I I t 1=1 t NO V ARS I 
20047  FORHATI//I*  *»TUPLtS«* * *101 / t Ix • GF 20 .1011 
30140  CONTINUE 

C  •••«* 

00  1200  I=1*N0VARS 
INPUTIII=TUPtJESIP  ill 
1200  CONTINUE 

CALL  PREPRI INPUT* VALUE  1 
TUPLES! P >N0 V All I = VALUE 
C  CHECK  FOR  A  SHAtL  GRADIENT 
SUM=0. 

DO  1300  I=1*M0VARS 
SUN=SUMi'GRA0(P-liII**2 
1300  CONTINUE 

SUH=SUH** I .51 
«=P-1 

IF  CSUH  .LT.  FLATI  0=0-1 
IF  (SUM  .LT.  FLATI  GOTO  3100 
C  TO  EXIT  IF  NO  SIBNIFICANT  CHANGE  FROM  PREVIOUS  VALUES 
00  3G00  I=1*N0VARS 

IF  (ABSITUPLE3HP*Il-TUPLESIP-l*II I  .GE.  EPSLONI  GOTO  3700 
3600  CONTINUE 
0=0-1 
GOTO  3100 
3700  CONTINUE 

C  END  OF  POINT  COMPUTATION  LOOP 
3000  CONTINUE 
3100  P=Q«1 

00  1400  Irl.NOVARS 
INPUT(II=TUPL£SIP ill 
1400  CONTINUE 

DO  1500  I=liNaVARS 

CALL  PARTLINeVARS*I*INPUT*EPSLON»PARTYL*INCR.TUPLES(P*NOVALlll 

C  •«••• 

IF  (SNITCH  .NE.  -II  GOTO  30150 

URITEI0UTFIL*Q003ai  (INPUT U I  * J= 1*N0V ARS I  * P ART YL * INCR * 

1  TUPLESIPiNOVALl  I 

20030  FORMAT!//*'  ••LAST  INPUT  PARTI AL^* * *10 1 / *1 X*6F20 .101  I 
30150  continue 
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GRADIPf  II  =  P*R7yL 
1500  CONTINUE 

IF  IQ  .LT.  II  Q=1 


IF  ISUITCH  .HE.  -1)  eOTO  30130 
URITE(0UTFILt20050l  CUM 
20050  FORNATI//t*  ••SUM • »F 30.15  I 
30130  CONTINUE 

C  ••••• 

C  begin  OUTPUT  OF  COMPUTATIONS 

C  WRITE  FILE  FOR  VARIABLE  VARIATION  ROUTINE  IVARYVARI 
WRITEIFILEltia02ai  P  tNOVARS.EPSLCN 
10020  FORMATI I5iI3tF10.7 I 

MRITEIFlLEliia030l  F0RMT2 
10030  FORMAT! ABOI 

00  3200  1  =  1  iP 

WRITEIFILE1.F0RMT21  ITUPLES 1 1 . Jl t J  =  1 . NOV  ARC  I . TUPLES  1 1 ,NOV ALl I 
3200  CONTINUE 
C  OUTPUT  TO  PRINIER 

MRITEIOUTFILiIQOAOI 
lOQitO  FORMAT!///. 21X.»TUPLES*t/l 
DO  3300  I=1.P 

WRITE !00TfIL.FCRHT3»  ! TUPLES! I »J » .J=l .NO WAPS  I . TUPLES  1 1 .MOV ALl I 
3300  CONTINUE 

WRI TEIOUTFIL. 100501 
10050  F0RMAT!///.2W.*GRA0IENTS*»/» 

DO  3400  1  =  1. P 

NRITCI0UTFIL.E0RMT3I  1 GR AD!I .U I .U=l .NOV ARS I 
3400  CONTINUE 

WRI TEIOUTFIL. 10060) 

10060  F0RMATI///.21X.* INCREMENT S*./» 

DO  3500  1=1. G 

WRI TEIOUTFIL. F0RMT3  I  ! DELTA! I, Jl . J=l .NO VARSI 
3500  CONTINUE 

WRITE!FILE2.10070i  P  .NOVARS.EPSLON 
10070  F0RNATII5.I3.F10.7) 

MRITEIFILE2.10030I  F0RMT2 
DO  3800  1  =  1. P 

NRITEIFILE2.F0RMT2I  IGRAD I I. J) . J=l. NOVARS) 

3800  CONTINUE 

WRITEIOUTFIL.IQOSO) 

10090  FORMAT! ///I 
END 
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2-8.  POINTCOMP  ROUTINE  FLOWCHART 


Read  in  format 
Tor  disk  and  / 
format  for  i 
printed  outnuti 
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Section  II.  THE  VARVARYl  ROUTINE 


2-9.  INTRODUCTION.  The  purpose  of  this  routine  is  to  take  each 
of  the  points  generated  by  the  POINTCOMP  routine  or  chosen  by  the 
user  and  to  vary  each  variable  in  turn  (using  the  bound  and  incre¬ 
ment)  while  keeping  the  other  variables  fixed.  The  program  to  be 
tested  is  evaluated  at  each  of  these  new  points,  and  the  new 
points  and  their  associated  program  values  and  gradients  are 
printed  out.  Difference  quotients  are  computed  and  printed  for 
each  variable  as  it  is  varied.  The  points  at  which  the  difference 
quotient  is  computed  are  those  defined  by  the  user  provided  incre¬ 
ments  and  lower  bounds.  The  new  points  can  be  used  to  study  the 
effects  of  varying  just  one  variable  at  a  time  and  to  provide 
points  for  graphs. 

2-10.  LIMITATIONS 

a.  When  varying  a  single  variable,  the  only  values  tested  are 
those  of  the  form:  lower  bound  +  integer  x  increment  i  higher 
bound.  The  current  routine  will  not  accept  a  list  of  values  of 
the  form  1,  3,  5,  17,  98. 

b.  All  variables  are  real. 

c.  Currently  restricted  to  100  variables. 

d.  Currently  restricted  to  a  maximum  of  500  steps  per  variable 
being  varied. 

e.  When  this  routine  is  used  on  output  from  POINTCOMP  or  GRID, 
the  limitations  of  these  routines  apply  as  well. 

2-11.  RUN  SETUP 

a.  Developing  an  Absolute  File  in  ASCI 
0MAP,S  ,  name  of  absolute  element 

IN  03PR0GTEST. VARVARYl 
IN  03PR0GTEST. PARTIAL 
IN  element  containing  PREPR 
IN  program  to  be  tested 
LIB$*FTN.8. 
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b.  To  Execute 


@LISE  10,  name  of  file  containing  input  points 

@USE  11,  name  of  file  into  which  gradients  are  placed  in 
standard  format 

@USE  12,  name  of  file  into  which  new  points  are  placed  in 
standard  format 

@USE  13,  name  of  scratch  file 

(3USE  14,  name  of  scratch  file 

@ASG,A  name  of  file  containing  input  points 

@ASG,A  name  of  file  into  which  gradients  will  be  placed 

@ASG,A  name  of  file  into  which  new  points  are  placed 

@ASG,A  name  of  scratch  file 

@ASG,A  name  of  scratch  file 

@XQT  absolute  element 

[Input  deck] 

c.  Description  of  Input  Deck 

(1)  Line  1.  (Format  to  read  in  each  of  the  following  three 
lines,  one  at  a  time.) 

(2)  Line  2.  Lower  bounds  for  the  variables;  the  leftmost 
bound  pertains  to  the  first  variable.  The  line  must  adhere  to  the 
above  format. 

(3)  Line  3.  Upper  bounds  for  the  variables. 

(4)  Line  4.  Increments  for  the  variables. 

(5)  Line  5.  (The  format  to  print  out  the  point  identifica¬ 
tion  number,  the  new  variable  valuer,  and  the  output  value  from  the 
routine  being  tested,  one  point  at  a  time.)  This  format  is  also 
used  to  print  out  the  gradient  identification  number  and  the  gra¬ 
dient  coordinates. 
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(6)  Line  6.  (The  format  to  print  out  one  line  of  difference 
quotients."] 

d.  Sample  Input  Data  for  VARVARYl 

(13F5.2)  Format  for  reading  in  the  next 

three  lines,  one  at  a  time. 


5.0 

0.5 

0.5 

0.5 

Lower  bounds. 

18.0 

9.0 

48.0 

1.0 

Upper  bounds. 

3.5 

1.5 

10.0 

0.1 

Increments 

The  leftmost  column  in  rows  1-4  pertains  to  the  first  variable, 
the  middle  column  to  the  second  variable,  etc. 

e.  Sample  Input  File  for  VARVARYl— Unit  10 

('  POINT  NUMBER  ' , I5,2X,8F10.5)  Format  for  printing  out  new 

points  and  gradients,  one  to  a 
1  i  ne. 

(1X,13F10.5)  Format  to  print  out  one  line  of  difference  quo¬ 
tients.' 


INPUT  VALUES 


9  4  .0010000 

(6(7X,F13.5)) 


5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

.61902 

.75809 

.63333 

.99916 

5.00119 

.73820 

.97516 

.76667 

1.24159 

5.00270 

.80924 

1.08781 

.84667 

1.37902 

5.00270 

.84061 

1.13401 

.88222 

1.43833 

5.00270 

.85178 

1.14995 

.89492 

1.45926 

5.00270 

.85513 

1.15468 

.89873 

1.46551 

5.00270 

.85513 

1.15590 

.89873 

1.46626 

5.00270 

.85513 

1.15715 

.89975 

1.46753 

This  file  was  created  by  POINTCOMP  on  unit  10. 


(1)  Line  1.  The  first  number  is  the  number  of  points  in  the 
15  format  and  the  second  is  the  number  of  variables  in  an  13  for¬ 
mat.  The  third  number  is  used  to  detennine  when  a  number  is  es¬ 
sentially  zero,  in  F10.7  format. 
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(2)  Line  2.  The  format  with  which  the  following  lines  were 
written,  and  with  which  they  may  be  read. 

(3)  Lines  3-7.  The  leftmost  four  columns  represent  the  val¬ 
ues  of  variables  one  to  four,  reading  from  left  to  right.  Each 
entry  in  the  rightmost  column  gives  the  tested  program  value  when 
run  with  the  variable  values  printed  on  the  same  row.  Each  row 
represents  an  input  point.  The  leftmost  four  columns  are  values 
of  the  input  variable.  The  last  column  is  the  output  value  of  the 
routine  being  evaluated. 

f.  Sample  VARVARYl  Run  Setup 

@USE  10,03MAT1. 

@ASG,A  03MAT1. 

@USE  11,03MAT3. 

@ASG,A  03MAT3. 

@USE  12,03MAT4. 

@ASG,A  Q3MAT4. 

@USE  13,03MAT5. 

@ASG,A  03MAT5. 

@USE  14,03MAT6. 

@ASG,A  03MAT6. 

@XQT  03PR0GTEST. VARVARYl 

Input  cards; 

(13F5.2) 

5.0  0.5  0.5  0.5 

18.0  9.0  48.0  1.0 
3.5  1.5  10.0  0.1 

('  POINT  NUMBER  ' , I5,2X,8F10.5) 

(1X,13F10.5) 
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2-12.  OUTPUT  DESCRIPTION  AND  SAMPLE  OUTPUT 
a.  Printer 


(1)  Sample  VARVARYl  Output— Points 


VARYING  VARIABLE  NUMBER  1 


POINT  NUMBER  1 
POINT  NUMBER  2 
POINT  NUMBER  3 
POINT  NUMBER  4 
POINT  NUMBER  5 


5.00000  .50000 
8.50000  .50000 
12.00000  .50000 
15.50000  .50000 
19.00000  .50000 


.50000 

.50000 

.50000 

.50000 

.50000 


.50000 

.50000 

.50000 

.50000 

.50000 


.73821 

.73934 

.73985 

.74014 

.74033 


++++++++++++++++++++DIFFERENCE  QUOTIENT  MATRIX+++++++++++++++++++ 


.00000 

.00032 

.00023 

.00018 

.00015 


.00000 

.00015 

.00011 

.00009 


.00000 

.00008 

.00007 


.00000 

.00005  .00000 


VARYING  VARIABLE  NUMBER  2 


POINT  NUMBER 

6 

5.00000 

.50000 

.50000 

.50000 

.73821 

POINT  NUMBER 

7 

5.00000 

2.00000 

.50000 

.50000 

1.10039 

POINT  NUMBER 

8 

5.00000 

3.50000 

.50000 

.50000 

1.32998 

POINT  NUMBER 

9 

5.00000 

5.00000 

.50000 

.50000 

1.48843 

POINT  NUMBER 

10 

5.00000 

6.50000 

.50000 

.50000 

1.60434 

POINT  NUMBER 

11 

5.00000 

8.00000 

.50000 

.50000 

1.69280 

POINT  NUMBER 

12 

5.00000 

9.50000 

.50000 

.50000 

1.76252 

h++++++++++++++++++D I FFERENCE 

QUOT I  ENT  MATR I X+++++++++++++++++++ 

.00000 


.24145 

.00000 

.19726 

.15306 

.00000 

.16672 

.12935 

.10564 

.00000 

.00000 

.14436 

.11199 

.09145 

.07727 

.00000 

.12728 

.09874 

.08063 

.06812 

.05897 
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(2)  Sample  VARVARYl  Output— Gradients 


THE 

VARIABLE  VARIED 

IS  NUMBER 

1 

POINT 

NUMBER 

1 

.00048 

.30919 

.67048 

.34638 

POINT 

NUMBER 

2 

.00019 

.32471 

.68200 

.34773 

POINT 

NUMBER 

3 

.00010 

.33081 

.68706 

.34836 

POINT 

NUMBER 

4 

.00006 

.33495 

.68991 

.34871 

POINT 

NUMBER 

5 

.00004 

.33762 

.69173 

.34894 

THE 

VARIABLE  VARIED 

IS  NUMBER 

2 

POINT  NUMBER 

6 

.00048 

.30919 

.67048 

.34638 

POINT  NUMBER 

7 

.02158 

.18565 

.60644 

1.08123 

POINT  NUMBER 

8 

.05104 

.12313 

.56450 

1.55265 

POINT  NUMBER 

9 

.07778 

.08812 

.53491 

1.88027 

POINT  NUMBER 

10 

.10199 

.06543 

.51291 

2.12195 

POINT  NUMBER 

11 

.12182 

.05001 

.49591 

2.30725 

POINT  NUMBER 

12 

.13869 

.04014 

.48239 

2.45384 

(3)  Description  of  Points  Output.  Each  row  describes  a 
point.  The  leftmost  four  columns  are  values  of  variables  one 
through  four,  reading  from  left  to  right.  The  rightmost  entry  on 
each  row  is  the  output  value  when  the  inputs  are  those  in  columns 
one  through  four. 

(4)  Description  of  Gradients  Output.  Each  row  is  the  gradi¬ 
ent  of  the  routine  being  tested  at  the  point  whose  number  is 
listed  on  the  left.  Again,  the  values  of  variables  one  to  four 
are  listed  from  left  to  right.  The  point  numbers  link  the  print¬ 
out  of  the  points  and  variable  values  to  the  gradient  printout. 

(5)  Sample  VARVARYl  Output—Difference  Quotients 


GIVEN  SAMPLE  OUTPUT 


POINT  NUMBER  51 

5.00000 

.73820 

.97516 

.76667 

1.24158 

POINT  NUMBER  52 

8.50000 

.73820 

.97516 

.76667 

1.26175 

POINT  NUMBER  53 

12.00000 

.73820 

.97516 

.76667 

1.27087 

POINT  NUMBER  54 

15.50000 

.73820 

.97516 

.76667 

1.27607 

POINT  NUMBER  55 

19.00000 

.73820 

.97516 

.76667 

1.27942 
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We  define  the  following  labels  for  some  of  the  variable  values  and 
outputs : 


Values  Outputs 


X,  =  5.00000 

Xp  =  8.50000 

Xo  =  12.00000 

x:  =  15.50000 
X5  =  19.00000 


Y,  =  1.24158 
Yp  =  1.26175 
Y^  =  1.27087 
Y4  =  1.27607 
Y5  =  1.27942 


THE  ASSOCIATED  DIFFERENCE  QUOTIENTS  ARE 


00000 

00576 

.00000 

00418 

.00261 

.00000 

00328 

.00204 

.00148 

.00000 

00270 

.00168 

.00122 

.00096 

.00000 

If  we  denote  the  varying  coordinate  of  point  numbers  51-55  by  x, 
to  xj,  and  denote  the  corresponding  values  (the  rightmost  column) 
by  yj  to  y^,  the  second  row  in  the  matrix  has  the  nonzero  entry: 

H  -  >1 

-  =  .00576. 

X2  -  xi 

Note  that  the  difference  in  y  is  in  the  numerator  and  the  differ¬ 
ence  in  X  is  in  the  denominator.  The  third  row  has  the  nonzero 
entries: 


^3  -  yi 

-  =  .00418 

^3  -  M 


^3  ■  H 

-  =  .00261 

X3  -  Xp 

Note  that  the  third  row  has  the  first  term  of  numerators  and  de¬ 
nominators  indexed  by  3  and  the  second  terms  are  indexed  by  1  and 
2  (i.e.,  3-1).  The  fourth  row  has  the  nonzero  entries: 
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y4  -  yi 

-  -  H  =  .00328 

X4  -  xi 


^4  -  H 

X4  -  X2 


.00204 


^4  -  ^3 

X4  -  X3 


.00148 


Note  that  the  fourth  row  has  the  first  term  in  each  numerator  and 
denominator  indexed  by  4  and  the  second  terms  are  indexed  by  1,  2, 
and  3  (i.e.,  4-1).  In  general,  the  nth  row  is  comprised  of  the 
ordered  set: 


3"*^  *. .  #  ,n— 1 ) 


To  show  the’  geometrical  significance 
the  following  points: 

of  these  computations 

given 

POINT  NUMBER  51 

5.00000 

.73820 

.97516 

.76667 

1.24158 

POINT  NUMBER  52 

8.50000 

.73820 

.97516 

.76667 

1.26175 

POINT  NUMBER  53 

12.00000 

.73820 

.97516 

.76667 

1.27087 

POINT  NUMBER  54 

12.50000 

.73820 

.97516 

.76667 

1.27607 

POINT  NUMBER  55 

19.00000 

.73820 

.97516 

.76667 
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and  the  difference  quotients: 


.00328 

.00270 


.00204 

.00168 


aflflOOO 

LooTTir 

1.00122 


Table  2-1.  Table  of  Difference  Quotients 


The  purpose  of  these  computations  will  now  be  explained.  The  im¬ 
port  of  the  quotients  in  triangle  A  is  illustrated  in  Figure  2-1 
(drawn  not  to  scale). 


slope  .00261 

slope  .00576  / 

\  ^  ^  slope  .00418 


53  54 


Figure  2-1.  Triangle  A  in  Table  2-1 
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The  figures  in  triangle  B  show  that; 


slope  .00148 


I 


slope  .00096 

slope  .00122 


_ u 


53 


54 


55 


56 


Figure  2-2.  Triangle  B  in  Table  2-1 


The  fourth  line  gives  the  following  picture  (not  to  scale); 


slope  .00204- 


slope  .00148 


A 

slope  .00328 


51 


52 


53 


- 

54 


„ .  -f  - 

55 


Figure  2-3.  Slopes  on  the  Fourth  Line  of  Table  2-1 
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Note  that  since  this  is  the  fourth  row,  all  lines  have  the  fourth 
point  as  terminus,  reading  from  left  to  right.  The  second  column 
gives  the  following  picture  (not  to  scale): 


slope  .00261 


slope  .00204 


slope  .00168 


51 


— 

52 


53 


1 

54 


55 


Figure  2-4.  Slopes  on  the  Second  Column  of  Table  2-1 

Note  that  since  this  is  the  second  column,  all  lines  start  at 
point  number  52,  reading  from  left  to  right. 
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225  4  .0010000 


(6(7X,F13.5)) 

5.00000 

.50000 

.50000 

.50000 

.73821 

8.50000 

.50000 

.50000 

.50000 

.73934 

12.00000 

.50000 

.50000 

.50000 

.73985 

15.50000 

.50000 

.50000 

.50000 

.74014 

19.00000 

.50000 

.50000 

.50000 

.74033 

The  first  row  indicates  that  there  are  225  points  in  the  file, 
that  there  are  four  variables  (the  rightmost  column  gives  the  out¬ 
put  values),  and  that  numbers  smaller  than  .001  are  considered  to 
be  zero.  The  second  row  gives  the  format  in  which  the  file  was 
written,  which  may  be  used  for  reading  in  the  file.  The  225 
points  and  values  commence  at  line  3  and  comprise  the  remainder  of 
the  file.  (Only  five  points  are  illustrated.) 

(2)  Standard  Format  File  on  Unit  11— Gradients 


225  4  .0010000 

(6(7X.F13.5)) 


.00048 

.30919 

.67048 

.34638 

.00019 

.32471 

.68200 

.34773 

.00010 

.33081 

.68706 

.34836 

.00006 

.33495 

.68991 

.34871 

.00004 

.33762 

.69173 

.34894 

The  first  two  rows 
dients  comprise  the 

of  this  file  are 
following  rows. 

described  as 

above.  The 

(3)  Scratch  Files 
(a)  Unit  13 


-19.00100 

4.00100 

.00048 

.30919 

.67048 

.34638 

20.00100 

4.00100 

.00055 

.37265 

.66454 

.34684 

21.00100 

4.00100 

.00061 

.43536 

.65859 

.34731 

22.00100 

4.00100 

.00067 

.49823 

.65264 

.34778 

23.00100 

4.00100 

.00072 

.56128 

.64668 

.34825 

24.00100 

4.00100 

.00076 

.62451 

.64071 

.34872 

25.00100 

4.00100 

.00079 

.68^90 

.63473 

.34920 

-26.00100 

1.00100 

.00359 

.35940 

.65458 

.40207 

CAA-D-80-1 

The  first  column  is  a  point  count.  A  negative  sign  indicates  a 
new  variable  is  being  varied.  The  second  column  indicates  which 
variable  is  being  varied.  Columns  3  through  6  are  the  gradient 
values  for  variables  1  to  4,  respectively. 

(b)  Unit  14 


5.00000 

.50000 

.50000 

.50000 

.73821 

8.50000 

.50000 

.50000 

.50000 

.73934 

12.00000 

.50000 

.50000 

.50000 

.73985 

15.50000 

.50000 

.50000 

.50000 

.74014 

19.00000 

.50000 

.50000 

.50000 

.74033 

5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

2.00000 

.50000 

.50000 

1.10039 

first  four 

columns  indicate 

variable 

values.  The 

fifth  column 

gives  the  output  values,  each  one  corresponding  to  the  input  val¬ 
ues  listed  on  the  same  row. 
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2-13.  VARVARYl  ROUTINE  LISTING 


P»RAHETER  F IL E 1= 10 . ^ IL E4 -1 3 . NO V AL S-IOD t INFILE = 5 . OUTFILrG 
PARAKETER  POINTErSQO  .F ILE 2  =  1 1 . F rL£3r 12 , F ILE 5 rlV 
C 

DIMENSION  TUPLE  S<  NOV  ALE J . LC WE OS « NOV AL S » . HYBOS ( NOV ALS I . INCR ( NO V ALS) 
DIMENSION  VALUrSENC V4LSI » GR AD « NO V AL S I  » V ALU f PC  INTS  »  »L I NE I  POINT S J 
C 

real  TUPLES.LOWBOS.HYBDStlNCR.VALUr S. VALUF. VALU.LINE.INCRl 
real  PART YL fF StFO  .  tPSLCN 
C 

INTEGER  NOV ARS tNOPIS  .1 t J .K  .L  .P tC 
C 

CHARACTFR*90  F C RM T 1 . F Q fi NT  2 . FORMTS  .  F0PMT4 
C 

REWIND  FILEI 
REWIND  FILE2 
REWIND  files 
REWIND  FILEv 
REWIND  files 
0  =  0 

REA D< FILE IflDOOD  J  NOF IS . N OV ARS tFPSLON 
lUOUO  formate I5»I3tF10. 71 

READEFILEl.lliniOJ  fqrmTI 
lOOlO  FORMATEASO) 

REACE  INFIlC  flUGlOE  FC'>MT2 

READEINFlLF,roFMT21  ELOWBOSEII t 1=1 . NO VARSI 
PEACE  infile  tFCPMT2  »  E H YSDS E I ) »I= ItNOV ARS ) 

REAOE INFIlE .F0RMT21  E INCR E 1 », 1=1 t NO VARS » 

READE  INFIlF  .liitllUl  F0RMT3 
REAOE  iNFlLf  »lunil3)  F0RMT4 
WRI TE  EOUTFiLilOnSUI 
luoao  FORMATE /////) 

DO  4U0  P=1.N0PTS 

REAOE file ItFCPM Til  E T UFL E S E I J , 1= 1,N0V*RS» 

CO  3LIU  I=liNOVAPS 
WRI TE EOUTFIL il0U70 »  I 

1D07D  FCRMATElXtllitlH-l . 'VARYING  VARTABLF  NUMBER  • . 1 3 . 10 E IM- > t // I 
DO  500  J=liNOVARS 
VALUESE JI=TUPLESEJI 
500  CONTINUE 

100  K  =  INTI I HYBDSI I )-L0W«C5Er  I  l/INCREI  |  1  +  1 

CO  2U0  J=OtK 

V  ALUESE  U=L0WBDSE  !»♦  J»TNCREI » 

CALL  PRFPRE V»LUFS fVALUE  I 
0  =  0'^! 

WRITE  EOUTFIL »F0RM Til  G . E V ALUE  S EL  I t L  =  1 .NO VA RS I .VALUE 
WRITE E files .FCPMTi  1  E V ALUES E L I .L  =  1 .NOVARS) .VALUE 
VALUE J*1I=VALUF 
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FQ-Q+EPSLON 
rX= I+EPSLON 
DO  150  l;i.nov«rs 

CALL  P ART L«  NO VARS •L . VALUES .EP SLOW t P APT YL  f INCR 1 » VALUE  1 
GRAOtD^PAKTYL 
15D  CONTINUE 

IF  (J  .Eg.  D)  EQr-r.l 

WRITECrlLE4»F0RMTl»  FO  .F I  ,  « GR AD J L > .L -1 , NO V A RS I 
20U  continue 

WRITE lOUTFiL. 100501 

10050  FORMAT!/, IX. 2QHH+Jt‘0IFFERENCE  QUOTIENT  M  ATR IX  •  ,  20  ( 1 H*)  ,/ ) 
00  700  J=0,K 
00  600  L:0,J 

LINE!  L  +  H  =  !  V  ALU!  J  +  ll  -VALU!L*1  »  »  / !  H  A  XO !  J- L ,  1 1  •  INCR  1 1>  » 
cm  CONTINUE 

WRITE  !OUTFILi  FORM  T4  1  (  LINE  IL  1 1 ,  L-U  .  J  > 

700  CONTINUE 

WRITE!0UTFlL«ia06DJ 
IDOGO  FORMAT! /( 

3U0  CONTINUE 

WRITE IOOTF1L.1D030>  P 

10030  FORMAT! /,  IX  ,  ?ll !  IH/ », ‘END  OF  COMPUTATIONS  FOR  POINT  NUMBER  •» 
1  I5.3X.30!1H/) .///> 

400  CONTINUE 

WRITE  !  FILE!  ,10000  »  0  ,NOV  ARS  ,  EP'^LON 

WRITE  !FlLE2,in010l  "^CRMTi 

WRI TE ! files ,10000  I  Q  tNCVAflS.EPSLCN 

WRITE IFILES ,100107  FORMTl 

REWIND  FILE4 

REWIND  files 

WRITE!OUTFIL, 100901 

10070  FORMAT!//, IX, 3U!1H  I, ’THE  CORRESPONDING  GRADIENTS*,/! 

DO  SOU  P:l,0 

READ! FI LE 4. FORM ri  I  Fa.FI,l5R&DILI,L:l,N0vARS! 

PEAO! file  5, FORM Til  ! V « LUE 3 ! L ) ,1: 1 , NOV A RS ) , V AL UE 
IF  IFG  ,GT.  0.1  GOTO  750 
F  Or-FQ 
K=INT!FTi 

WRITEIOUTFlL.iniOOl  K 

10100  FORMAT! /, 10! IH  I,*  THE  VARIABLE  VARIED  IS  NUMBER  *,15,/) 

750  J=INT!FQ) 

WRITE ! OUT FIL, FORM T3  I  J , ! 6R AC I L I , L  =  1 , NO V A RS  » 

WRITE!FlLE3, FORMTl!  !  V ALuE S I L I ,L:l , NO V AR S!  .VALUE 
WRITE ! file 2, FORM! II  !GR AO ! L  I  ,L r 1 , NO V ARS 1 
900  CONTINUE 

WRITEIOUTFIL, 100401 
10040  FORMAT!///! 

END 
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CHAPTER  3 

THE  ANALYZ  SUBSYSTEM 


3-1.  INTRODUCTION.  This  routine  computes  variable  sensitivity 
statistics  from  gradients  in  a  standard  format  file  on  unit  11. 
Statistics  are  initially  computed  on  each  gradient.  The  gradients 
are  summed,  component  by  component,  and  the  same  statistics  are 
computed  for  the  sum  vector. 

3-2.  DISCUSSION  OF  STATISTICS 

a.  The  first  statistic  is  the  ratio  of  each  component 

of  the  gradient  to  the  component  with  the  smallest  abso¬ 
lute  value. 

b.  The  second  statistic  is  the  change  in  the  component  re¬ 
quired  to  achieve  a  unit  change  in  the  output  for  each  component 
of  the  gradient.  For  the  summed  vector,  a  change  in  the  output 
equal  to  the  number  of  gradients  rather  than  a  unit  change  is 
utilized. 

c.  The  third  statistic  is  the  ratio  of  each  component  to  the 
component  with  the  largest  absolute  value;  except  that  for  the 
largest  component,  its  ratio  with  the  next  largest  component  is 
taken. 


3-3.  LIMITATIONS 

a.  The  gradients  input  to  ANALYZ  must  all  be  real. 

b.  No  more  tiian  500  gradients  can  be  analyzed  in  a  single  run 
by  the  currently  compiled  version  of  ANALYZ. 

c.  No  more  than  20  variables  may  be  analyzed  by  the  current 
version  of  ANALYZ. 

3-4.  RUN  SETUPS 

a.  To  Execute 

yuSE  11,  name  of  file  containing  gradients  in  standard  format. 
@ASG,A  filename. 

@XQT  03PR0GTEST. ANALYZ 
[Input  deck] 
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b.  Description  of  Input  Deck 

(1)  Line  1.  Field;  A1 

Format:  13 

A1  is  +1  if  no  debugging  information  is  desired,  -1  if  debugging 
is  desired. 

{?)  Line  2.  (Format  for  printing  the  statistics.) 

(3)  Line  3.  (Format  for  debug  printouts.) 

c.  Sample  Input  Data 


+1  Indicates  no  debugging  output  re¬ 

quested. 

(6(7X,F13.5),/,100(21X,5(F13.5,7X),/))  Format  for  statisti¬ 

cal  output 

(6(7X,F13.5),/,100(21X,5(F13.5,7X),/))  Debugging  data  output 

format 


d.  Input  File  Description  on  Unit  11 


Produced  by  POINICOMP  on  Unit  11 

9  4  .0010000 

(P(7X,H3.b)) 


.00048 

.30919 

.67048 

.34638 

.00359 

.35940 

.65458 

.4020/ 

.00852 

.40086 

.63565 

.45144 

.01236 

.42295 

.62294 

.4/942 

.01428 

.43221 

.61698 

.49139 

.01500 

.43545 

.61481 

.49560 

.01522 

.43641 

.61415 

.49686 

.01523 

.43631 

.61415 

.49674 

.0152/ 

.43670 

.61405 

.49663 

(1)  Line  1 
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(a)  Format:  15,  13,  F10.7. 
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(b)  Meaning:  Nine  lines,  four  variables,  ^ero  approxima¬ 
tor  .001. 

{Z)  Line  2.  Format  for  reading  in  the  following  lines  of 
data,  one  at  a  time. 

(J)  Other  Lines.  Each  row  is  a  gradient,  variable  1  values 
are  in  the  leftmost  column,  variable  2  values  are  in  the  next  col¬ 
umn,  etc. 

e.  Sample  Run  Setup 

@USt  11,03MAT2. 

@ASG,A  03MAT2. 

@XQT  03PR0GTEST.ANALYZ 
+1 

(6(7X.H3.b),/,100(21X,5(F13.5,7X),/)) 

(6(7X,F13.5),/,100{21X,5(F13.5,7X),/)) 

3-b.  OUTPUT  DESCRIPTION  AND  SAMPLE  OUTPUT 
a.  Sample  Output--Gradient  bj/  Gradient 

STATISTICS  ON  EACH  GRADIENT 

***********’*******«*'****[:gy  j  Q^^^Q^^************************ 


.00000 

3.23426 

278.55153 

2.78242 

117.37089 

2.49464 

80.90615 

2.36435 

70.02801 

2.31369 

66.66667 

2.29647 

65.70302 

2.29142 

65.65988 

2.29195 

65.48788 

2.28990 

1.48147 

2.88/00 

1.52770 

2.48713 

1.57319 

2.21513 

1.60529 

2.08585 

1.62080 

2.03504 

1.62652 

2.01776 

1.62827 

2.01264 

1.62827 

2.01313 

1.62853 

2.01357 
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J  !!•****♦★★*■**’*★  !!«*★*★★★**★★★*★♦ 


.00000 

.46155 

1.93568 

.51661 

.00548 

.54905 

1.62802 

.61424 

.01340 

.63063 

1.40805 

.71020 

.01984 

.67896 

1.29936 

.76961 

.02314 

.70053 

1.25558 

.79644 

.02440 

.70827 

1.24054 

.80610 

.02478 

.71059 

1.23606 

.80902 

.02480 

.71043 

1.23636 

.80883 

.02487 

.71118 

1.23643 

.80878 

***it'kitieir'kiriciti(*icitir*4tirieir-k*ief^£[^l!^Jiy£  J  J  Qg*********************** 


.00000 

1.00000 

2.16850 

1.12028 

1.00000 

100.11142 

182.33426 

111.99721 

1.00000 

47.04930 

74.60681 

52.98591 

1.00000 

34.21926 

50.39968 

38.78803 

1.00000 

30.26681 

43.20588 

34.41106 

1.00000 

29.03000 

40.99733 

33.04000 

1.00000 

28.67346 

40.35151 

32.64520 

1.00000 

28.64806 

40.32502 

32.61589 

1.00000 

28.59856 

40.21284 

32.52325 

b.  Description.  The  preceding  statistics  are  described  in 
paragraph  3-?.  ^r  each  statistic,  each  row  corresponds  to  a 
gradient. 

c.  Sample  Output — Sum  of  Gradients 


SUMMLU  GRADIENT  COMPUTATIONS 


*★★***-**  J  ******************  ****** 

.09947  3.66948  5.65779  4.15653 


********************fQ^jjy/\LENT  SUMMED  CHANGES********************* 
90.47954  P.45266  1.59073  2.16427 
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***ir4r****'***************3mif||V|[;Q  QQf/jp^p  [  3Q(^5************************ 

.017b8  .64857  1.36118  .73466 


**********************pf SUMMED  RATIOS********************** 
l.OOOOn  36.89032  56.87936  41.78677 


d.  Description.  The  first  row  is  the  sum  of  the  gradients. 

The  last  three  rows  are  the  same  statistics  as  before,  in  the  same 
order,  but  applied  to  the  sum  vector  only. 


3-5 


l,AA-D-80-l 


36.  ANALYZ  ROUT.  -  LISTING 


1. 

?. 

3. 

4. 

5. 
S. 
7. 
3. 
3. 

lU. 
11. 
12. 
1  3. 
14. 
13. 
16. 
17. 

19. 
13. 

20. 
21. 
7  2. 

23. 

24. 
?S. 
26. 
27. 

29. 
23. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 
3  7. 
39. 
3  3. 
4U. 


PARAMETEP  FILE2rll,:NETLr-5.CUTFIL=6 
parameter  N0VAL5:?U  .NOGRA'J- 5lin 

r 

PTMEN3I0N  GRAD  eNCGRAO.NOVALS  )  •RMINS  I  NOGPAD  J  . C '^uMS  ( NO V AL S  » 
or  KEN  51  ON  RKAX5(NOGRAO,2).CG.MAX5(?>  .PRTNT  I  NOV  AL  S  7 
C 

P^AL  grad .RMING .C50M5. "KAXS  »C5MAX5*PR  ENT. EP3L0N 
C 

INTEGER  C  .I.V  .f.‘0VAR3.N0t.P!>.  IWITCM 
C 

CHARACTEP *8U  FKTRO.FmTWRT.FMTCE 1 

c 

C  READ  the  inputs 
C 

REAOI  infile  t  llhl)  SWITGH 
iCiU  FORMAT!  13  1 

R^AD! INFTLE .iUUI  FMTU3T 
2Ulj  FORMAKASO) 

C 

C  READ  IN  THE  GRA0?''Nt  VALU'D 
C 

REWIND  FILE2 

READ! FI LE2. SOU  I  NC5 RO . NO V AR S . EF SLO N 
3UIj  forma  T!  IE  ,13.F10.7  1 

READ! FILE2. 200)  FMTRD 
DO  4110  GrliNCGRO 

READ! FILE2.FMTR0)  !GFA0!G.V).V:l.N0VA35) 

400  continue 

C 

C  READ  IN  dcbug  mode  FORMAT  IF  IN  rFRU”^  MODE 
C 

IE  (SWITCH  .EG.  -1)  ^EADCINFILE ,?DUI  FMTD3G 
C 

IF  (SWITCH  .NE.  -1)  GOTO  20200 
WRITE (OUTFiLtZUOPOl  NOCPD.NOVAFS 
20(I0U  FORMAT!//.  •  •* I NPUT  ) 

00  2LltlU  C:l.NOGPD 

WRITE  (OUTFIL.FM’GPG  I  (  ■"  R  AD  (  G  .  V  )  .  1 .  NO  V  AR  S  ) 

/niOlI  CONTINUE 

20. Of)  continue 


4  1. 
4  2. 
4  3. 
4  4. 
45. 
«6. 
4  7. 
48. 
43. 


C  . 

r 

C  COMPUTE  ROW  MINIMA 
C 

DO  750  G-l.NCGRC 
CO  500  V=l.NrV*RS 
RMINS(GI=0.0 

IF  (A9S (GPAD(C.VI)  .IT.  FPSLON)  GOTO  SOU 
RMINS(G l  =  AgS( GRAD(G .Vt  ) 
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5lJ.  GOTO  oUU 

f.i.  soil  continue 

52.  tOU  00  7U0  V:liNOV*RS 

53.  IF  IftBStCRADtGiV >1  .UT.  EPSLON)  PC AC  I G . V)  =  U. 0 

54.  IF  (AR3(CRA0(G.Vn  .LT.  EPSLOM»  GCTC  TOO 

55.  IF  ( ABS IGRAOl G| V >  I  .UT.  RMlNSCCit  RNINS ( C ) =AB3 I GRAO ( G  •  V  I  I 

5G.  700  continue 

57.  750  CONTINUE 

58.  C  ••*»«.»••• 

59.  TF  (SWITCH  .NE.  -H  COTC  20400 

ED.  WRITE  l0UTFlL«2U30ai 

61.  203011  FOPMAK//.’  ••RCW  MlUS**‘t/> 

C2.  WRITE  (OUTFiLtFMTOP.GI  (  RMINS  I  G  J  i  C  =  1 »  NOf  R  D I 

6  3.  2UI.U11  CONTINUE 

fA.  C 

F  5.  C 

66.  C  COMPUTE  THE  COLUMN  SUMS 
E7.  C 

6  3.  00  3U0  V:l*NOV<fRS 

69.  CSUMStVJrfi. 

7D.  CO  aOU  Grl.NOGRO 

71.  CSUMS(V»:CSUKS(WI+G«AD(G,V( 

72.  SQL  CONTINUE 

73.  IF  <AaS<CSUMS(Vn  .LT.  IN0G«0»EPSU0N)  )  CSUMStVl=D.D 

74.  90U  continue 

75.  CSMN:AaS(C3UHSIll  I 

76.  CO  lUOU  V^Z^NCVARC 

77.  IF  (4aS<CSUHS ( VI  I  .LT.  ( NCCRO*r P5L0N n  GOTO  IDUD 

78.  IF  (ABStCSUMStVU  .LT.  CSMINJ  C3MIN=  A33  K  SUM9  W  I  1 

79.  1000  CONTINUE 

80.  C 

81.  C  TC  COMPUTE  THE  TWC  OICTINCT  LARGEST  FARTIALSIIN  ARSOLUTE  VALUE! 

92.  C  IN  EACH  GRADIENT 

8  3.  C 

84.  IF  (NOVARS  .LT.  21  GOTi  L6UU 

85.  no  1350  G-l.NCGRD 

86.  PMAXS(G,1 |:ABS(GRA0 (C ,11  I 

97.  DO  IIUU  V:2.N0VARS 

89.  IF  (  ABS  I  GPAO(  G,  V  )  I  .CT.  RMAXSIC.D)  RM  A  XS  t  G 1 1 )  :  A8  S  (  GRAD  t  6  .  V  11 

89.  iiOU  CONTINUE 

90.  DO  IISU  V^I.NCVARC 

91.  IF  |ABS(GRA0(G.V1  I  .LT.  EP5L0NJ  GOTO  1150 

32.  IF  (A03(r,RA0(G,vn  .CE,  RMAXSIG,!))  GOTO  1150 

93.  RMAXS(Gt2l-AeS(GR4D(G.VI » 

94.  GOTO  1200 

95.  1150  CONTINUE 

96.  RHAXS(G«2 I^RMAXS (Of  1  I 

97.  GOTO  1350 
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93.  120D  PO  13UO  V=1.KCV&93 

99.  IF  (  »  AE  3«  RRAf.' (G  .V  J  »  .f  T  .  R«  AX  3  I  G  t  ?  )  »  .AND.  » A93  •  GRAD  «  G  .  V  I  »  .LT. 

IDO.  1  RMAXSIGtDll  R“  A  X  S  C  f  A93  «  GR  An  (  G  t  V  »  ) 

ini.  1300  CONTINUE 

102.  1350  CONTINUE 

103.  C  »*»..*«.«* 

lOM.  IE  ICWITCH  .NE.  -1!  CGTo  -jogU 

IDS.  WRITE  IOUTEIL«21l7nni 

ins.  21)700  FORMATI/.*  *»ROU  MAXIMA. «♦,/) 

107.  DO  2030U  G=1.N0CR0 

10  9.  WRITE  •OUTFIL.EMTDFG  I  «  r  I?  AP  <  G  .  V  I  .  1 .  NO  V  AR  C  )  »  PM  A  XS  <  G  .  1 »  .RMAXS  «  G.  2» 

109.  209U0  CONTINUE 

110.  21000  CONTINUE 

111.  C  .».»..»»*. 

112.  C 

113.  C  TO  CCMFUTE  THE  TWO  LARGEST  'THPONENTS  OF  CSUHS  IN  ARSOLUTE  VALUE 

im.  c 

115.  n^iMAXSI  1  I-A95  ( '^UMS  (  U  1 

IIG.  DO  1300  V-2tN0VARS 

111.  IF  (  A9'J  (CGUMC  (  V  I  J  .07.  C-^MAXSti))  C  CM  A  X  S  1 1  )  =  A9S  (  C  SU  MS  «  V  I  I 

J13.  1300  CONTINUE 

119.  no  1325  Vrl.NOVARS 

120.  IF  ( A GS CCSUMS< V  I  1  .LT.  C N^C RP »E P SL ON )  I  GOTO  1325 

121.  IF  ( APS (CSUMS 7 V t  I  .  GF  .  CSMAxS<ll)  GOTO  142  5 

122.  rSMAXS(2>rA9SIC  ,UM3«V)  » 

123.  GOTO  1350 

12  3.  13^5  continue 

125.  C^MAXSlPIrCSMAX^dl 

12S.  COiO  1550 

127.  1450  no  1500  Vrl.NCVARS 

123.  If  (  (  AB'H  CSUHSI  V  M  .''T.  CSMAXSID)  .AND.  <  ABE  ( C3UMS  I  V  )  I  .LT. 

121.  I  rSMAXS(l)ll  C  SMAXS<2  »r  APS«CSLM''(  V  )  » 

130.  150U  CONTINUE 

131.  155U  CONTINUE 

132.  C 

133.  IF  C-HITCM  .nE.  -II  rjTO  21300 

133.  WRITE (0UTFIL.2110UI 

135.  211UU  FORMAT!/.*  ••CSUMS  MAyI9A.»*,/l 

13S.  WPTTE(CUTFlL«FMTn9GI  IC':U!'S!VI,V:l»N0VARS),CPMAxSm.nSMAXSI2» 

197.  215UU  continue 

133.  C 

139.  IbUU  CONTINUE 

130.  C 

131.  C  TO  COMPUTE  THE  RELATIVE  PATIC: 

13  2.  C 

131.  WRITE  lOUTEIL.lt'OOOl 

133.  lOUOU  FORMAT!//, IX. 3L(1H»I,*PFLATIV  R A T I  OS  ’  . ID  I 1H» ) , / I 
135.  PO  leou  Grl.NOoRC 
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me. 

m?. 

1<|9. 

149. 

150. 

151. 

152. 
15  1. 
164. 
155. 
15£. 
157. 
15  8. 
159. 
lEO. 
IE  1. 
IE  2. 
IE  3. 
ie4. 
ie5. 
lEE. 
1E7. 
If  8. 
If  3. 

170. 

171. 

172. 

17  3. 

174. 

175. 

176. 

177. 

178. 

179. 

180. 
181. 
182. 

18  3. 
184. 
IBS. 
166. 
187. 
183. 

189. 

190. 

191. 

192. 


00  17CJII  V  =  l.NOVAi?S 
PRINTIV)=n.0 

IF  I  RMIN^CGI  .3F.  EP5LON  I  PR INH V  I :GR AD t G t V ) /RMTNS ( G I 
1700  CONTINUE 

WRITE  COUTFILtFHTWRT  I  « PP TNT < V I . V: 1 • NO VARS  I 
1800  CONTINUE 
C 

C  PRINT  THE  equivalent  CHANC'' 

C 

WRITE  «OUTFlLtlOinC» 

lolun  F0fiHATI//tlX»3miH*l. •EQUIVALENT  CH  ft  NGF  S  •  t  3u  C IH  »  I  ,/ » 

DO  2000  C^ltNOCfiO 
DO  13U0  V^l.NCVARS 
PRINT (V  1  =  0.0 

IF  1  GRADIG.V)  .GF.  rPSLON  )  PR  IN T ( V ) = 1 . /G RAP ( 6 . V ) 

1900  CONTINUE 

WRT  TE  «0UTFIL*FMTWRT  1  « F PINT • V »  1 7=1 . NO VARS ) 

2000  CONTINUE 
C 

C  PRINT  THE  COMPARISONS 
C 

IF  (NCVARS  .LT.  21  3CT0  2250 
WRI TE (OUTFIL. 102001 

10201!  rOPMA  r(//.  IX.  IUHH»  }  F  •COMPARISONS’  .ID  I IH*  I  ,  /  » 

00  2200  Prl.NOGRO 
DO  21U0  V=1.N0VARS 
PRINT(V»:0.Q 

IF  (  RMAXS<3.1I  .GE.  EPSOON  I  PRINTI VJrGRADI G.V)/flMAXS(G»l ) 

IF  ( A6S JPRINT  •  V  )  )  .EG.  !•>  PRINT ( V ) :GR AO  I G  .  V ) /RMAXS C G .2 » 

210U  CONTINUE 

WRITE  (OUTFIL.FMTWPT  1  (  PR  ttlT  <  V  »  .  V  =  1 .  NO  VARS  » 

2200  CONTINUE 
2250  CONTINUE 
C 

C  SUMMED  ORADirNT  CCMPUTATTON. 

C 

WRITE(CUTFIL.lli3nnj 

linuo  PORMATI//.IX. 31 (IH+ I .’CUMMEO  GRAOTENT  COMPUT ATI0N3' ,301 1H+ I ./) 
C 

C  WRITE  OUT  THE  SUMMED  GRADIENTS 
C 

WRITFI0UTFIL.1035U) 

10350  rOPMATI//. IX, 3li(lH«l .’SUMMED  OR APTE NT S ’ , 30  1 1H» J , / I 
WRITE (OUTFIL.FMTWPT •  ( C SUMS • V »  ,  V: 1 . NO  VARS  I 
C 

C  COMPUTE  THE  SUMMED  RATIOS 
C 
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1S3.  WRITF  IOUTFIL*104tin» 

194.  lU40tl  FORM*T» // »1X.  3U«  IH*  »  .  •RFUTiVE  SU*'MFD  P  AT lOS  *  »  3U  (  iH*  I  . /  » 

195.  DO  lUSDU  V=1.NCVARS 

19G.  P91NT(V):n.U 

197.  IF  (  CSMIN  .GE.  rpsLON  »  PR INT I  V  I : C SUMS « V t /CGMI N 

193.  lUf.DU  CONTINUE 

199.  WRI  TE  (OUTFIUFMTWRT  I  ( '^N  UV  )  .  V  -  1 .  NO  V  ARS  » 

?UO.  C 

201.  C  COf’PUTr  THE  SUMMED  EGUI  VALIN!  CHANCES 

202.  C 

203.  WRI TE (OUTFIL.lOfcnOl 

204.  lUEUQ  <^ORMAT( //,1X.  Till  IH.  I  ,  JL'’'VALENI  SUMMED  C  H  AN  -  ES  *  .  30  « IH*  I  »/ ) 

205.  DO  lOTOU  V-l.NCVAPS 

2P6.  FRINTCVim.U 

207.  IF  (  CSUM'‘(VI  .GE.  EP^LON  I  PRl  NT  I  V  )  rl.'OGPD/CGUMS  ( V  I 

20  3.  1I)7PL'  CONTINUE 

209.  WRI Tr (oUTFILtFMTWRT »  1 p  IlNT » V > » V T 1  . NO V ARG J 

21U.  C 

211.  C  COMPUTE  SUMMER  COME  ARISONS 

21  2,  C 

213.  IF  (  NOVARS  .LT.  21  !‘070  11000 

214.  WRI TE (OUTEIl. lOaOO) 

215.  luauU  •'ORMATI //.IX.  31. « 1H»  ).  •SUMMED  COMP  AR  IS  0  N  S  '  .  30  ( 1 H  *  I  .  /  » 

216.  DO  lUatiO  V:l.NtVAPS 

217.  PRINHVJrO.O 

213.  IF  (  CSMAXSdl  .GE.  EPSLC.N  )  PR  I N  T  (  V  >  :  C  GUMS «  V  )  /  CSM  A  XS  U  » 

219.  TE  (ABS(PRINT« VI  I  .EO.  1.1  PRI N T ( V  I ; C SUMS t V I / CSM A XS C 2  I 

220.  10900  CONTINUE 

221.  WRITE IOUTFIL.fhTWRT  I  < PS INT I  V 1 . Vr 1 , NO  VARS  I 

222.  linOO  CONTINUE 

.’23.  WPI  Ti- (OUTFIL.llltjlll 

224.  llinu  ECRMATI///I 

225.  ENC 
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CHAPTER  4 

THE  GRID  SUBSYSTEM 


Section  I.  THE  GRID  ROUTINE 

4-1.  INTRODUCTION.  This  routine  utilizes  a  user  specified  sub¬ 
division  of  each  variable  in  order  to  generate  a  grid  of  input 
variable  values.  The  routine  to  be  tested  is  evaluated  at  each 
node  of  the  grid,  and  the  gradients  are  also  computed  at  each 
node.  The  output  of  the  grid  routine  is  used  by  REARRANGE  and 
DIFFQUOT. 

4-2.  LIMITATIONS 

a.  As  cu,  rently  compiled,  the  routine  is  restricted  to  20  in¬ 
put  variables  at  the  most. 

b.  All  input  variables  (being  varied)  must  be  real. 

c.  Only  one  output  from  the  routine  being  tested  may  be 
checked  out  at  one  time. 

d.  Testing  time-consuming  Monte  Carlo  routines  may  be  too  ex¬ 
pensive.  Rather  than  making  several  runs  per  point  and  averaging 
the  output  values,  it  is  better  to  run  one  component  of  the  system 
at  a  time,  testing  random  variables  no  differently  from  determi¬ 
nistic  variables. 

e.  The  function  represented  by  the  routine  being  tested  must 
be  well-behaved. 

4-3.  RUN  SETUPS 

To  Develop  an  Absolute  ASCI  Program  File 
@MAP,S  name  to  be  given  to  absolute  element 
IN  03PR0GTEST.GRID. 

IN  03PR0GTEST. PARTIAL. 

IN  element  containing  the  driver  PREPR. 

IN  programs  to  be  tested. 

LIB$*FTN8. 

END. 
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b.  To  Execute 


0USE  10,  name  of  file  into  which  points  will  be  stored. 

@ASG,A  name  of  this  file. 

0USE  11,  name  of  file  into  which  gradients  will  be  stored. 

PASG,A  name  of  this  file. 

OUSE  12,  name  of  scratch  file  (used  as  input  to  REARRANGE  and 

DIFFQUOT). 

0ASG,A  name  of  this  file. 

OUSE  13,  name  of  scratch  file. 

0ASG,A  name  of  this  file. 

OXQT  name  of  absolute  deck  created  by  OMAP. 

[Input  deck] 

c .  Description  of  the  Input  Deck 

(1)  Line  1.  Number  of  variables  -  13 

Zero  approximatot  -  F10.7 
Debug  mode  field  -  13 

The  zero  approximator  is  a  threshold-,  numbers  smaller  in  absolute 
value  are  considered  to  be  zero.  If  the  debug  mode  field  contains 
the  numhet  -1,  the  run  will  be  in  debug  mode.  Any  other  value  im¬ 
plies  run  will  not  print  debugging  infoniiation. 

(2)  Line  2.  (Format  for  reading  in  one  line  of  variable 
variat  ion  tTa ta. ) 

(3)  Line  3.  Initial  variable  values.  Variable  one's 
initial  value  is  leftmost,  variable  two’s  initial  value  is  to  the 
right  of  variable  one,  etc. 

(4)  Line  4.  Bounds  for  variable  values. 

(5)  Line  5.  Increment  (step)  value  for  each  variable. 
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(6)  Line  6.  Fonnat  for  one  line  of  output.  First  field 
should  be  16,  other  fields  should  be  real.  The  number  of  fields 
specified  must  be  no  less  than  number  of  variables  +  1. 

(7)  Line/.  File  I/O  fonnat  for  one  line  of  output.  All 
fields  should  "he  real  and  at  least  (number  of  variables  +  2) 
fields  must  be  specified. 

(8)  Line  8.  Optional;  should  contain  a  fonnat  for  writing 
one  line  of  debuggi ng  data  if  in  debug  mode. 


d.  Sample  Input  Data 
4  .001  +1 


(13F6.0) 


Four  variables.  Zero  ap¬ 
proximator  value  of  .001. 
Not  in  debug  mode 


Format  for  reading  vari¬ 
able  variation  data 


5.0  0.5  0.5  0.5  Initial  values 


18,00  9.00  48.00  0.9 


3.6  1.5  10.0  O.J 


Terminating  values  for  in¬ 
crementation 

Step  values  for  each  vari¬ 
able 


('  POINT  NUMBER  M6.5(7X,F13.51))  Format  for  output 


(6(7X,F13.5)) 


File  I/O  format 


Not  a  debugging  run,  so  debugging  format  is  omitted, 
e.  Sample  Run  Setup 
WUSE  10,03MAT5. 

OASG.A  03MAT5. 

OUSE  11,03MAT6. 

OASG.A  03MAT6. 

OUSE  12,03MAT7. 

OASG.A  03MAr7. 
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@USE  13,03MAT8. 

@ASG,A  03MAT8. 

@XQT  03PR0GTEST.C0MGRID1. 


4 

(13F6.0) 

.001 

+1 

5.0 

0.5 

0.5 

0.5 

18.00 

9.00 

48.00 

0.9 

3.5 

1.5 

10.0 

0.1 

('  POINT  NUMBER  '  ,I6,5(7X,F13.5)) 

(6(7X,F13.5)) 

4-4.  OUTPUT  DESCRIPTIONS  AND  SAMPLE  OUTPUT.  These  outputs  are  of 
the  same  types  as  the  outputs  of  POINTCOMP  and  are  comprised  of 
variable  values  and  gradients.  The  two  routines  differ  in  the 
point  selection  algorithm. 

a.  Sample  Output--Variable  Values  and  Output  Values 
VARIABLE  NUMBER  3  HAS  BEEN  INCREMENTED 


POINT  NUMBER 

31 

5.00000 

.50000 

50.50000 

.50000 

34.26218 

POINT  NUMBER 

32 

5.00000 

.50000 

50.50000 

.60000 

33.99982 

POINT  NUMBER 

33 

5.00000 

.50000 

50.50000 

.70000 

33.73710 

POINT  NUMBER 

34 

5.00000 

.50000 

50.50000 

.80000 

33.47402 

PUINT  NUMBER 

35 

5.00000 

.50000 

50.50000 

.90000 

33.21060 

POINT  NUMBER 

36 

5.00000 

.50000 

50.50000 

1.00000 

32.94681 

VARIABLE  NUMBER  2 

HAS  BEEN 

INCREMENTED 

POINT  NUMBER 

3/ 

5.00000 

2.00000 

.50000 

.50000 

1.10039 

POINT  NUMBER 

38 

5.00000 

2.00000 

.50000 

.60000 

1.20860 

POINT  NUMBER 

39 

5.00000 

2.00000 

.50000 

.70000 

1.31727 

POINT  NUMBER 

40 

5.00000 

2.00000 

.50000 

.80000 

1.42641 

POINT  NUMBER 

41 

5.00000 

2.00000 

. 50000 

.90000 

1.53502 

POINT  NUMBER 

42 

5.00000 

2.00000 

.50000 

1.00000 

1.64610 
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VARIABLE  NUMBER  3  HAS  BEEN  INCREMENTED 


POINT 

NUMBER 

43 

5.00000 

2.00000 

10.50000 

.50000 

7.16476 

POINT 

NUMBER 

44 

5.00000 

2.00000 

10.50000 

.60000 

7.08334 

POINT 

NUMBER 

45 

5.00000 

2.00000 

10.50000 

.70000 

7.00158 

POINT 

NUMBER 

46 

5.00000 

2.00000 

10.50000 

.80000 

6.91946 

POINT 

NUMBER 

47 

5.00000 

2.00000 

10.50000 

.90000 

6.83699 

POINT 

NUMBER 

48 

5.00000 

2.00000 

10.50000 

1.00000 

6.75416 

POINT 

NUMBER 

48 

5.00000 

2.00000 

10.50000 

1.00000 

6.75416 

In  each  row,  the  four  columns  to  the  right  of  the  point  number 
contain  the  four  input  variable  values,  and  the  fifth  contains  the 
corresponding  output  value  from  tlie  routine  being  tested.  The 
first  variable,  as  usual,  is  the  leftmost,  i.e.,  in  the  column  to 
the  right  of  the  point  number.  Note  that  in  each  group,  the  last 
(fourth)  variable  is  being  varied  through  its  range  while  the 
other  variable  values  remain  fixed.  The  headings  explain  how  this 
group  differs  from  the  preceding  group.  Whenever  the  first,  sec¬ 
ond,  or  third  variable  is  incremented,  the  variables  to  its  right 
are  reset  to  their  initial  values. 


b.  Sample  Qutput--Gradients 


VARIABLE  NUMBER  3 

HAS  BEEN  INCREMENTED 

POINT 

NUMBER 

31 

.26871  -2.38501 

.67048 

-2.62300 

POINT 

NUMBER 

32 

.32286  -2.86548 

.66454 

-2.62653 

POINT 

NUMBER 

33 

.37827  -3.34741 

.65860 

-2.63007 

POINT 

NUMBER 

34 

.43280  -3.83228 

.65264 

-2.63362 

POINT 

NUMBER 

35 

.48479  -4.31738 

.64668 

-2.63718 

POINT 

NUMBER 

36 

.54242  -4.80347 

.64071 

-2.64073 

VARIABLE  NUMBER  2 

HAS  BEEN  INCREMENTED 

POINT 

NUMBER 

37 

.02158  .18565 

.60644 

1.08123 

POINT 

NUMBER 

38 

.02576  .22479 

.58747 

1.08586 

POINT 

NUMBER 

39 

.03090  .26340 

.56843 

1.09051 

POINT 

NUMBER 

40 

.03513  .30235 

.54930 

1.09520 

POINT 

NUMBER 

41 

.03931  .34164 

.53010 

1.09991 

POINT 

NUMBER 

42 

.04344  .38127 

.51081 

1.10466 

VARIABLE  NUMBER  3 

HAS  BEEN  INCREMENTED 

POINT 

NUMBER 

43 

.15683  -.14934 

.60644 

-.81352 

POINT 

NUMBER 

44 

.18874  -.17997 

.58747 

-.81700 

POINT 

NUMBER 

45 

.22084  -.21085 

.56843 

-.82051 

POINT 

NUMBER 

46 

.25426  -.24200 

.54930 

-.82403 

POINT 

NUMBER 

47 

.28688  -.27446 

.53010 

-.82758 

POINT 

NUMBER 

48 

.31969  -.30626 

.51081 

-.83115 
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The  point  numbers  link  the  gradients  to  the  variable  values  at 
which  the  gradients  were  computed. 

c.  Sample  Output  Files.  On  unit  10--points  in  standard  for¬ 
mat. 

1260  4  .001000 


(6(7X,F13.5)) 

5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

.50000 

.50000 

.60000 

.77276 

5.00000 

.50000 

.50000 

.70000 

.80735 

5.00000 

.50000 

.50000 

.80000 

.84199 

5.00000 

.50000 

.L.0D00 

.90000 

.87667 

5.00000 

.50000 

.50000 

1.00000 

.91140 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

The  first  line  is  the  standard  heading,  component  of: 

(1)  The  number  of  points  in  15  format. 

(2)  The  number  of  variables  in  13  format. 

(3)  The  zero  approximator  in  F10.7  format. 

The  second  line  is  the  format  used  to  write/read  the  file.  The 
following  lines  comprise  the  variable  values  and  corresponding 
output  values. 

d.  On  Unit  ll--Gradients  in  Standard  Format 


1260  4  .0010000 

(6(7X,F13.5)) 


.00048 

.30919 

.67048 

.34638 

.00055 

.37265 

.66454 

.34684 

.01,1061 

.43536 

.65859 

.34731 

.0006/ 

.49823 

.65264 

.34778 

.00072 

.56128 

.64668 

.34825 

.00076 

.62541 

.64071 

.34872 

.05325 

-.22734 

.67048 

-.24905 

.06396 

-.27318 

.66454 

-.24938 

.07469 

-.31915 

.65859 

-.24972 

The  first  two 

lines  are  as  described 

previously. 

The  following 

lines  contain 

gradient  values  corresponding  to  the 

variable  values 

contained  in 

the  previously  described 

file. 
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e.  On 

Unit  12--File 

Used  t  o 

Communicate 

with  DIFFQUOT  and 

REARRANUr 

5.00000 

.50000 

.50000 

.50000 

.73821 

4.00100 

5.00000 

.50000 

.50000 

.60000 

.77276 

4.00100 

5.00000 

.50000 

.50000 

.70000 

.80735 

4.00100 

5.00000 

.50000 

.50000 

.80000 

.84199 

4.00100 

5.00000 

.50000 

.50000 

.90000 

.87667 

4.00100 

5.00000 

.50000 

.50000 

1.00000 

.91140 

4.00100 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

3.00100 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

4.00100 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

4.00100 

In  each  row,  the  first  four  columns  represent  variable  values. 

The  fifth  column  contains  output  values.  The  integer  portion  of 
the  sixth  number  is  the  variable  being  varied  in  obtaining  the 
values  for  that  row.  The  rightmost  column  is  used  by  REARRANGE 
and  DIFFQUOT.  This  file  is  written  and  may  be  read  by  the  format 
located  at  the  second  line  of  the  standard  format  files  on  units 
10  and  11. 

f.  On  Unit  13--Scratch  File  Containing  Gradients 


00048 

.30919 

.67048 

.34638 

.73821 

00055 

.37265 

.66454 

.34684 

.77276 

00061 

.43536 

.65859 

.34731 

.80735 

00067 

.49823 

.65264 

.34778 

.84199 

00072 

.56128 

.64668 

.34825 

.87667 

00076 

.62451 

.64071 

.34872 

.91140 

05325 

-.22734 

.67048 

-.24905 

7.44301 

06396 

-.27318 

.66454 

-.24938 

7.41817 

07469 

-.31915 

.65859 

-.24972 

7.39330 

The  first  four  columns  contain  the  gradient  components.  The  fifth 
column  contains  the  variable  values  corresponding  to  the  input 
variable  values  at  which  the  gradient  was  evaluated.  This  file 
was  also  written  and  may  be  read  using  the  format  on  the  second 
line  of  the  files  on  units  10  and  11. 


CAA-l)-8U-l 


4-0.  GRID  RUUliNL  LISTING 


PARAMETER  I NE TLEr 5 . OUTriLrE . F ILE l=lDt FIL r2 r 11 1 F IL E3r 12 , FILE<#  =  13 
PARAMETER  N 0 V ALS: 20 . NO VALl rNO V AL S*1 * NO V A L2 -NO VAL1*1 
C 

DIMENSION  FSTVAL  CNCVALSI tLSTVALINOVALS) . STEPS  C NOV ALS  J . V ALU C NO VALSI 
DIMENSION  COUNTS  I  NOV ALS I .LINE! NOV ALT  I t OR  AD (NOV ALS I 

C 

INTEGER  NOV ARSt SWITCH, POINTS. COUNTS 
INTEGER  PTR.PTRl.V.P 
C 

real  FSTVAL.LSTVAL.STEPS.V ALU.?'’ ad,  value  .LINE 
REAL  EPSLON.PARTYL.INCR.FPTR 

c 

CHARACTER .80  FMTRO .FMTPRT .FMTFIL . FMT08G 
C 

C  INTTIAL12ATI0N 

r> 

REWIND  FILE! 

REWIND  FILE2 
REWIND  files 
REWIND  ‘’ILE4 

REAO( infile .10000  I  NO V ARS . EPSLON , SW I T CH 
10000  FORMAT! I3.riO,7.I3» 

REAO( infile .10010 (  FMTRO 
loom  FORMAT(AeO» 

REAO( INFILE .FMTRO  I  ( F ST VAL I  V J . V- 1. NO V ARS  > 

READ! INFIlE .FmTRD  »  ( L STVAL I  V J , V: 1, NOV ARS » 

READ( INFIlE .FMTRO  I  ( STEPS C V ». V:1 . NC VARS ) 

P0INTS=0 

P,EAO(  INFILE  .lOOlU)  FMTPRT 
READ! INFILE  .10010  I  FMTFIL 

IF  (SWITCH  .EQ.  -II  REACdN'-ILE.lOUlO)  FMTDRG 
DO  100  V:l.NOVARS 
VALU( VlrFSTVAL(VI 
100  CONTINUE 

C  *.«•*••••• 

IF  (SWITCH  .N£.  -II  CO’O  20000 
WRITEIOUTFIL. 100201 

10020  FORMAT! /.IT  ,2li(lH-).  *£000  P  R I  M  RU  T  S  • .  20  (  IH- I  ,  /  I 
WRITEIOUTFIL. 100251  NOVARS.FF  , LON. SWITCH 
10025  F ORMATI 1X,I3,F10.T,T3I 
WR1TE(PUTF1l.1003UI 

10030  FORMAT! /, IX . •»•  FIRST  VALUES  •••! 

WRITE (CUTFIL.FMTDBC I  ( Fst V AL ( V  I . V=1 . NO V A RS I 
WRITE  (OUTFIL.IOOAL'I 

10040  FORMAT! /.IX  LAST  VALUES  »*’l 

WRI TE (OUTFIL.FMTOBGI  ( LST V AL ( V  I . V- 1 . NO V  A RS I 
WRI TE (OUTFIL.IOOSU  I 
10050  FORMAT! /, IX. INCREMENTS  .*•1 
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WRITf  <OUTFIL.FMTDaGt  ( ':TEP?  <  V  1 1  V=l*  NO  V»RS» 

20000  CONTINUE 

c  •••**•***• 

FPTR=FL0»TI NOVARSI^EPSLON 
GOTO  500 
C 

C  TO  CONFUTE  THE  NEXT  SET  OF  VARIABLE  VALUES 
C 

200  DO  TOO  PTR=NOVARS fl. -1 
FTR1=PTP 

FPTR^FLCATI PTRI l+EPSLON 

IF  <VALUCPTRI  .LT.  LSTVALCPTRJ )  GOTO  400 
VALU«PTR»:FSTVAL«PTRI 
300  CONTINUE 
GOTO  600 

4  00  VALU(PTR1 1  - VALU<  PTRi » ♦ STE PS « P TR 1 1 
500  CONTINUE 

POINTS-POINTS+1 

CALL  PREPR( V«LUt VALUE » 

WRITEIFILE3#FMTFIL  I  IVALUC  VI  .Vn.NOVAfiSI  tVALUEtFPTR 
00  550  V=liNCVARS 

call  PARTLI N0VARS«V. VALU.EFSLONfPARTYLtlNCRtWALUEl 
GRAOl VirPARTYL 

C  «**•••**** 

IF  (SWITCH  .NE.  -II  GOTO  20010 
WRITE(0UTFILtl006UI  POINTS. V 

10060  FORMAT!/. IX. ’»•  PARTIAL  SUBROUTINE  OUTPUT  FOR  POINT  •, 

X  16. •  ANO  variable  ‘.la.AH  »») 

WRITE (OUTFIL.FMTDBGI  P ART YL . VALUE . INC R 
2U010  CONTINUE 

C  •«««*••••* 

5Sn  CONTINUE 

WRI TE (FILE4 .FmTFIL I  ( GR AD ( V I . V =1 . NO V ARS I  .VALUE 
C  RETURN  TO  COMPUTE  THE  NEXT  SET  OF  VARIARLF  VALUES 
GOTO  200 

600  REWIND  files 
REWIND  FILE4 

WRITEIFILEl.lOOTUI  PO INTS .NOV ARS .EPSLON 
1007C  FCRMATI I5.I3.F10.7I 

WRITECFiLEl.iUniOl  FMTFIL 
WRI TE(FrLE2 .10070  I  POINTS.NOVAPS.EFSLON 
URITE(FILE2.1001UI  FMTFIL 
WRrTE(0UTFIL.10075l  POINTS 

10075  FORMAT!///. 31X. 'THE  ♦.15.’  POINTS  AND  OUTPUT  VALUES'. /I 

DO  aon  p=i.poiNTs 

READ!  FILE  3.  FMTFIL!  !  L  INF  !  W  I  .  V- 1 .  NOV  ARS  I  .  LIN'^ !  NO  V  AL 11  .  L INE  I  NOV  AL2  I 
PTR--INT(LINE!N0VAL2I  I 
IF  (PTR  .EG.  NOVARSI  GCTO  700 
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WRITE(0UTFIL»lUO8tjl  PTR 

lOniU  FORMAT!/. IX. ’VaRlARLE  NUMBER  ‘.IS.*  HAS  3FFN  IN CRFM ENTEO •  . 

1  / 1 

700  WRITEIOUTFIL.FMTPRTI  P . ( L INE ( V ) . Vrl . NO V A RS I . L INE I  NOV ALU 
WRITEfFTLEl.FMTFIL)  «L INE I  V  I . V"! . NOVARS ) .LINE (NOV ALl  I 
600  CONTINUE 

WHITE (OUTFIL.lOnaSJ 

1U015  FORMAT!///. 31X. 'THE  CORRESPONDING  GRADIENTS*./) 

REWINC  FILET 
DO  30U  P::l. POINTS 

READ! FTLE3.FMTFIL  I  ( L INF ( V ) . Vr I . NO V A R3 ) . LINE C NO V ALl )  . L INE ( NO V AL2  ) 
PTR-TNT(LINE(NCVAL2I  ) 

IF  (PTR  .EG.  NOVARSI  GOTO  850 
URITE(0UTFIL.HIQ3UI  PTR 
85U  PEAOl FILER. FMTFIL  )  ( L INE ( V ) » V- 1 . NO V  APS ) 

WRITE (OUTFIL.FMTPRT 1  P . ( L INE ( V ) . V=1 . NO V A RS ) 

WR1TE(FILE2.FMTFILI  I L IME ( V  I . Vr 1 , NO VARS  I 
9U0  CONTINUE 

WRITE(0UTFIL.1U03UI 
10U30  FORMAT! ///) 

ENO 
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graiJient 

heading 


Section  II.  THE  REARRANGE  ROUTINE 


4-7.  INTRODUCTION.  The  output  of  the  GRID  routine  is  so  organ¬ 
ized  that  within  each  group  of  output  lines,  the  lines  agree  in 
all  variable  values  except  the  last,  i.e.,  the  last  variable  is 
varied  most  frequently.  The  REARRANGE  routine  reorders  the  out¬ 
put,  grouping  on  each  variable  in  turn,  except  the  last.  This  re¬ 
ordering  facilitates  the  analysis  of  the  output. 

4-8.  LIMITATIONS 

a.  The  currently  compiled  version  is  limited  to  20  variables. 

b.  Since  this  routine  utilizes  output  from  GRID,  the  limita¬ 
tions  of  the  GRID  routine  are  applicable. 

c.  This  routine  uses  much  I/O,  so  the  smaller  the  number  of 
points  produced  by  GRID,  the  better. 

4-9.  RUN  SETUPS 

a.  To  Execute 


OUSE  10,  name  of  standard  format  file  produced  by  GRID  on 
unit  10. 

@ASG,A  name  of  above  file. 

(3USE  12,  name  of  communication  file  produced  by  GRID  on  unit 

12. 


(?ASG,A  name  of  above  file. 

OUSE  14,  name  of  file  to  be  rearranged.  This  file  should  be 
the  nonstandard  format  file  produced  by  GRID  on  units  12 
(points)  or  13  (gradients). 

(3ASG,A  name  of  this  file. 

OUSE  15,  name  of  scratch  file  into  which  REARRANGE  places 
outputs. 


0ASG,A  name  of  scratch  file. 
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@XQT  name  of  absolute  file  developed  earlier. 

[Input  deck] 

b.  Description  and  Sample  of  Input  Deck 

(1)  Line  1.  ('  POINT  NUMBER  M5,2X,7(F10. 5, 5X),/21X,8(F10. 5, 5X)) 


This  format  wil!  bring  out  the  point  number  and  regrouped  data. 

(2)  Line  2.  (15,2X,7(F10.b ,5X) ./.21X,8(F10.5,5X) ) 

This  format  will  be  used  to  printout  the  file  of  reordered  data, 
c.  Sample  Input  File 
(1)  Unit.  10 
1260  4  .0010000 


(6(7X,F13.5)) 

5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

.50000 

.50000 

.60000 

.77276 

5.00000 

.50000 

.50000 

.70000 

.80735 

5.00000 

.50000 

.50000 

.80000 

.84199 

5.00000 

,50000 

.50000 

.90000 

.87667 

5.00000 

.50000 

.50000 

1.00000 

.91140 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

5.00000 

.50000 

10.50000 

.80000 

7.36839 

5.00000 

.50000 

10.50000 

.90000 

7.34346 

This  is  a  standard  format  file  of  variable  values  and  an  output 
value.  This  file  is  developed  by  GRID  on  unit  10,  and  is  de¬ 
scribed  in  the  GRID  documentation. 
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(2)  Unit  12 


5.00000 

.50000 

.50000 

.50000 

.73821 

4.00100 

5.00000 

.50000 

.50000 

.60000 

.77276 

4.00100 

5.00000 

.50000 

.50000 

.70000 

.80735 

4.00100 

5.00000 

.50000 

.50000 

.80000 

.84199 

4.00100 

5.00000 

.50000 

.50000 

.90000 

.87667 

4.00100 

5.00000 

.50000 

.50000 

1.00000 

.91140 

4.00100 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

3.00100 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

4.00100 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

4.00100 

5.00000 

.50000 

10.50000 

.80000 

7.36839 

4.00100 

5.00000 

.50000 

10.50000 

.90000 

7.34346 

4.00100 

This  communication  file  is  developed  by  GRID  on  unit  12  and  is  de¬ 
scribed  in  the  GRID  documentation. 


(3)  Unit  14.  This  is  the  file  to  be  rearranged.  It  may  be 
the  GRID  file  produced  on  unit  12  described  earlier  or  the  analo¬ 
gous  nonstandard  format  gradient  file  produced  by  GRID  on  unit  13. 
A  portion  of  the  gradient  file  follows; 


.00048 

.30919 

.67048 

.34638 

.73821 

,00055 

.37265 

.66454 

.34684 

.77276 

.00061 

.43536 

.65859 

.34731 

.80735 

.00067 

.49823 

.65264 

.34778 

.84199 

.00072 

.56128 

.64668 

.34825 

.87667 

.00076 

.62451 

.64071 

.34872 

.91140 

.05325 

-.22734 

.67048 

-.24905 

7.44301 

.06396 

-.27318 

.66454 

-.24938 

7.41817 

d.  Sample  Run  Stream 


(aUSE  10,03MAT5. 
PUASG.A  03MAT5. 
OUSE  12,03MAT7. 
OASG.A  03MAT7. 
OUSE  14,03MAT7. 
OASG.A  03MAT7. 
OUSE  15,03MAT9. 
OASG.A  03MAT9. 


4-18 


@XQT  03PROGTEST. REARRANGE. 

(■  POINT  NUMBER  ' , I5,2X,7(F10.5.5X) ,/,21X ,8(F10.5 ,5X) ) 
(I5,2X,7(F10.5,5X),1.21X,8(F10.5,5X)). 


4-10.  OUTPUT  DESCRIPTIONS  AND  SAMPLE  OUTPUT 


a.  Printed  Output 


++++++1-+++++++++++++++++V  AR I ABL  E 

NUMBER 

1++++++++++++++++++++++++ 

POINT  NUMBER 

1 

.00048 

.30919 

.67048 

.34638 

.73821 

POINT  NUMBER 

253 

.00019 

.32471 

.68200 

.34773 

.73934 

POINT  NUMBER 

505 

.00010 

33081 

.68706 

.34836 

.73985 

POINT  NUMBER 

757 

.00006 

.33495 

.68991 

.34871 

.74014 

POINT  NUMBER 

1009 

.00004 

.33762 

.69173 

.34894 

.74033 

END  OF  VARIATION 

FOR  ] 

INITIAL  POINT  1 

POINT  NUMBER 

2 

.00055 

.37265 

.66454 

.34684 

.77276 

POINT  NUMBER 

254 

.00022 

.38998 

.67838 

.34892 

.77406 

POINT  NUMBER 

506 

.00012 

.39826 

.68446 

.34856 

.77465 

POINT  NUMBER 

758 

.00007 

.40214 

.68788 

.34887 

.77498 

POINT  NUMBER 

1010 

.00005 

.40531 

.69006 

.34908 

.77520 

END  OF  VARIATION  FOR  INITIAL  POINT  2 

This  particular  output  file  is  a  file  of  gradients,  a  reordered 
version  of  the  gradient  file  produced  by  GRID.  The  gradients  in 
the  first  group  were  evaluated  at  points  whose  values  for  vari¬ 
ables  2-4  are  identical.  Adjacent  gradients  were  evaluated  at 
points  whose  variable  1  coordinates  differ  by  the  step  value  for 
variable  1.  The  gradients  in  the  second  group  are  similar,  these 
gradients  were  also  evaluated  at  points  whose  coordinates  2 
through  4  are  identical  and  where  adjacent  gradients  were  evalu¬ 
ated  at  points  differing  only  in  the  variable  1  coordinate  and  the 
difference  is  the  step  size.  Note  that  the  fifth  column  contains 
the  corresponding  output  values.  The  points  at  which  these  gradi¬ 
ents  were  evaluated  are  exhibited  below. 
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++++++++++++++++++++++++VAR I ABLE  NUMBER  1++++++++++++++++++++++++ 


POINT  NUMBER 

1 

5.00000 

.50000 

.50000 

.50000 

.73821 

POINT  NUMBER 

253 

8.50000 

.50000 

.50000 

.50000 

.73934 

POINT  NUMBER 

505 

12.00000 

.50000 

.50000 

.50000 

.73985 

POINT  NUMBER 

757 

15.50000 

.50000 

.50000 

.50000 

.74014 

POINT  NUMBER 

1009 

19.00000 

.50000 

.50000 

.50000 

.74033 

END  OF  VARIATION 

FOR 

INITIAL  POINT  1 

POINT  NUMBER 

2 

5.00000 

.50000 

.50000 

.60000 

.77276 

POINT  NUMBER 

254 

8.50000 

.50000 

.50000 

.60000 

.77406 

POINT  NUMBER 

506 

12.00000 

.50000 

.50000 

.60000 

.77465 

POINT  NUMBER 

758 

15.50000 

.50000 

.50000 

.60000 

.77498 

POINT  NUMBER 

1010 

19.00000 

.50000 

.50000 

.60000 

.77520 

END  OF  VARIATION  FOR  INITIAL  POINT  2 


Note  that  within  each  group,  only  variable  1  varies.  The  follow¬ 
ing  printout  is  a  later  segment  of  the  same  output. 

FINISHED  VARIABLE  NUMBER  1 

++++++++++++++++++++++++VAR I ABLE  NUMBE R  2++++++++++++++++++++ ++++ 
NEW  INITIAL  POINT 


POINT  NUMBER 

1 

.00048 

.30819 

.67048 

.34638 

.73821 

POINT  NUMBER 

37 

.02158 

.18565 

.60644 

1.08123 

1.10039 

POINT  NUMBER 

73 

.05104 

.12313 

.56450 

1.55265 

1.32998 

POINT  NUMBER 

109 

.07^78 

.08812 

.53491 

1.88027 

1.48843 

POINT  NUMBER 

145 

.10199 

.06543 

.51291 

2.12195 

1.60434 

POINT  NUMBER 

181 

.12181 

.05001 

.49491 

2.30725 

1.69280 

POINT  NUMBER 

217 

.13869 

.04014 

.48239 

2.45384 

1.76252 

NEW  INITIAL  POINT 


POINT  NUMBER 

253 

.00019 

.32471 

.68200 

.34773 

.73934 

POINT  NUMBER 

289 

.00997 

.23380 

.63786 

1.18644 

1.15472 

POINT  NUMBER 

325 

.02581 

.17577 

.60447 

1.81117 

1.46142 

POINT  NUMBER 

361 

.04306 

.13756 

.57832 

2.29508 

1.69706 

POINT  NUMBER 

397 

.05990 

.10961 

.55728 

2.67928 

1.88374 

POINT  NUMBER 

433 

.07674 

.09007 

.54000 

2.99334 

2.03525 

POINT  NUMBER 

469 

.09131 

.07532 

.52555 

3.25328 

2.16066 

4-20 


CAA-D-80-1 


These  gradients  were  evaluated  on  point  groups  where  only  variable 
2  varied  within  each  group.  The  “NEW  INITIAL  POINT"  message  indi¬ 
cates  that  variable  1  varied  between  points  217  and  253,  so  these 
points  vary  in  two  coordinates— hence  it's  time  to  start  a  new 
grouping.  This  situation  will  not  occur  for  variable  1  because 
variable  1  is  varied  the  least.  The  following  printout  exhibits 
the  points  at  which  the  gradients  were  evaluated. 

FINISHED  VARIABLE  NUMBER  1 

++++++++++++++++++++++++VAR I ABLE  NUMBER  2++++++++++++++++++++++++ 


NEW  INITIAL  POINT 


POINT 

NUMBER 

1 

5.00000 

.50000 

.50000 

.50000 

.73821 

POINT 

NUMBER 

37 

5.00000 

2.00000 

.50000 

.50000 

1.10037 

POINT 

NUMBER 

73 

5.00000 

3.50000 

.50000 

.50000 

1.32998 

POINT 

NUMBER 

109 

5.00000 

5.00000 

.50000 

.50000 

1.48843 

POINT 

NUMBER 

145 

5.00000 

6.50000 

.50000 

.50000 

1,60434 

POINT 

NUMBER 

181 

5.00000 

8.00000 

.50000 

.50000 

1.69280 

POINT 

NUMBER 

217 

5.00000 

9.50000 

.50000 

.50000 

1.76252 

NEW  INITIAL  POINT 


POINT 

NUMBER 

253 

8.50000 

.50000 

.50000 

.50000 

.73934 

POINT 

NUMBER 

289 

8.50000 

2.00000 

.50000 

.50000 

1.15472 

POINT 

NUMBER 

325 

8.50000 

3.50000 

.50000 

.50000 

1.46142 

POINT 

NUMBER 

361 

8.50000 

5.00000 

.50000 

.50000 

1.69706 

POINT 

NUMBER 

397 

8,50000 

6.50000 

.50000 

.50000 

1.88374 

POINT 

NUMBER 

433 

8.50000 

8.00000 

.50000 

.50000 

2.03525 

POINT 

NUMBER 

469 

8.50000 

9.50000 

.50000 

.50000 

2.16066 

Note  that  all  points  in  each  group  vary  only  in  the  second 
variable.  Note  also  that,  as  before,  while  the  point  numbers  con¬ 
tinue  increasing,  points  217  and  253  differ  in  two  coordinates 
since  variables  1  and  2  change  simultaneously.  This  fact  necessi¬ 
tates  creating  a  new  group. 

b.  Sample  Output  File  Segment 


252 

5.00000 

9.50000 

50.50000 

1.00000 

15.70270 

504 

8.50000 

9.50000 

50.50000 

1.00000 

21.01859 

756 

12.00000 

9.50000 

50.50000 

1.00000 

24.55868 

1008 

15.50000 

9.50000 

50.50000 

1.00000 

27.08553 

1260 

19.00000 

9.50000 

50.50000 

1.00000 

28.97974 

252+++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
1  //////////////////////////////////////////////////// 
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VARIABLE  2 


NEW  INITIAL  POINT 


1 

5.00000 

.50000 

.50000 

.50000 

.73821 

37 

5.00000 

2.00000 

.50000 

.50000 

1.10039 

73 

5.00000 

3.50000 

.50000 

.50000 

1.32998 

109 

5.00000 

5.00000 

.50000 

.50000 

1.48843 

145 

5.00000 

6.50000 

.50000 

.50000 

1.60434 

181 

5.00000 

8.00000 

.50000 

.50000 

1.69280 

217 

5.00000 

9.50000 

.50000 

.50000 

1.76252 

The  leftmost  integers  are  point  numbers. 

( 1 )  The  1 i ne 

252  +  +  +  .  .  .  indicates  the  end  of  the  grouping  whose 
initial  point  was  252,  i.e.,  if  there  was  a  point  253, 
the  point  numbers,  now  at  1260,  would  decrease  to  253 
next  and  commence  to  increase  from  that  value,  i.e.,  the 
next  point  numbers  would  be; 

253 
289 
325 
etc. 

(2)  The  line 

1  ////  .  .  .  indicates  the  end  of  regroupings  for  which 

only  variable  1  varies  within  each  preceding  group. 

(3)  The  line 

VARIABLE  2  indicates  that  now  groupings  where  only 
variable  2  varies  will  be  derived. 

(4)  The  line 

- .  .  . - new  initial  point . indicates 

that  while  the  point  numbers  may  continue  to  increase,  a 
new  grouping  must  nevertheless  begin. 
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4-11.  REARkANGE  ROUTINE  LISTING 


F Afi AMETER  I NF ILE -  5 . OUT F I L  =  6 . F ILF ir 10 » FIL E3 =1 2 . SRT F IL i 14 . GRPFIL  =  15 
PARAMFfR  NOVAL':  =  20.NQVAL1:MOVALS*1*NOVAL2-NOVALS*2 
C 

CHARACTER»8U  FM TF IL t FMTPfi T .FHTGRP 
C 

INTEGER  POTNTS.PTRtSTEP.STFPN.PFRIOO.F.I.Q ,Q1,P,V 
INTEGER  NOVRMl tNO VARSfNOVRPl tN0VPP2 
C 

REAL  LINE 
C 

OIHENSTON  period ( NOV ALS»*LINE(NCVAL2 I 
C 

C  INITIALIZATION 
C 

REWIND  FILEI 
REWIND  FILE3 
REWIND  SRTFIL 
REWIND  GRPFIt 

READ(FILEl.lOQOQ)  POINTS  t  NO  VAR>; 
lOUOO  FORMAT! 15.13) 

REAOIFlLEl.lPOin)  FHTFIU 
loom  FORMAT!  AeO) 

READ! INFIlE .10010)  FMIPRT 
READ! iNFiLE .lOniO  )  FMTGRP 
N0VRP1=NCVARE+1 
N0VRP2=N0VARS+2 
N0VRM1:N0VAPS-1 

c 

C  TO  FIND  THF  PERIODS 
C 

DO  IbO  V:i.NOVSRS 
PEPICDC  V)-U 
100  CONTINUE 

DO  200  P;i. POINTS 

READ! FILE3.FMTFIL)  IL INE I  I » . 1  =  I . N0VRP2 ) 

PTR=INTILINEIN0VRP2I ) 

IF!  PERICDCPTR")  .EQ.  Ql  PERIOD!  PTR  I  rMAXDI  P- 1 . 1 ) 

IFlPERIOOd)  .NE.  PI  GOTO  300 
200  CONTINUE 
C 

C  TO  PRODUCE  A  REORDERED  FILE.  VARYING  EACH  VARIABLE  THRCUCHOUT 
r  TTS  range  in  TURN 

c 

300  DO  600  F^I.NOVRMI 

WRITE lOUTFiL. 100201  F 

10020  FORMAT! //.IX. 30(1H*I  .'VARIABLE  NUMBER  » , 15 . 2X  .  30  1 1 )  , // ) 

WRITE  GRPFlL.inOIU)  F 
10030  FORMAT! • variable  '.15) 
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STEPN=0 

STEP=PE(»ICn«FJ 

IF  IF  .NE.  II  STEPN:PERIOOIF-il 
C 

C  TO  PICK  UP  THE  Q  TH  ENTRY  FROH  EACH  BLOCK  OF  SIZE  STEP 
C 

DO  SOD  0=liSTEP 
REWIND  SRTFIL 
G1=Q 

IF  IQl  .EQ.  STEPI  ai=0 
CO  400  P=1»P0INTS 

PEAOl  SRTFIL  tFMTFlLI  ILINE  I V  I  *  V::l  t  NOVRPl  I 
IF  IMOOIP«STEP)  .NE.  Oil  GOTO  400 
IF  IQl  .EQ.  Ill  GOTO  350 
C 

C  CHFCK  FOR  VARIATION  IN  THE  NEXT  VAPTABLE 
C 

IF  IIF  .NE.  II  .AND.  INOOlPtSTEPNI  .EO.  0111  WRITE  lOOTFILtlOOJZI 
10032  rORMATI/.*  NEW  INITIAL  POINT*. /I 

IF  IIF  .NE.  II  .AND.  IHOOIP.STEPNI  .EQ.  0111  WRITE IGRPFIL . 10036 1 
10036  F0fiMATIS5llH-1.*NEW  INITIAL  POINT  * . 55  I IH-I  I 
350  WRITFIOUTFXL.FPTPRTI  P.IL1N£IVI.V=1»N0VARSI»LINEIN0WRP1I 
WRITEIGRPFIL.FKTGRPI  P » I  LINE  I  V I » V=1 .NOV ARS I .LINE  I NCVRPII 
IF  I  G1  .NE.  0  I  GOTO  400 

IF  IIF  .NE.  II  .AND.  IMODIP.STEPNI  .EQ.  Oil  WRITE lOUTFIL. 100321 
IF  IIF  .NE.  II  .AND.  I  MOO  IP .STEPN I  .EQ.  Oil  WRITE  I GRPFIL . 10036 1 
400  CONTINUE 

WPITEIOUTFIL. 100401  3 

10040  FORMATI/.*  ENC  OF  VARIATION  FOR  INITIAL  POINT  *.15.^1 
WRITEIGRPFIL. 100501  Q 
10050  FORMATI T5.120I1H»I  I 
C 

C  FINISHED  PASS  FOR  THE  0  TH  ENTRY  IN  EACH  BLOCK 
C 

SOU  CONTINUE 

WRITE lOUTFiL. 100601  F 

10060  FORMATI/.*  FINISHED  VARIABLE  NUMBER  *.131 
WRITEIGPPFILtlOOTOI  F 
10U70  F0RMATII5.5X.120I1H/II 
5UU  CONTINUE 

WRITEIOUTFIL.1003UI 
10030  FORMATI///I 
END 
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Section  III.  THE  DIFFQUOT  ROUTINE 

4-13.  INTRODUCTION.  This  routine  utilizes  the  output  developed 
by  GRID  on  unit  10  in  order  to  compute  difference  quotients  for 
one  variable.  The  point  numbers  listed  in  the  output  identify  the 
points  used  to  compute  the  difference  quotients.  These  points  are 
developed  by  GRID,  and  the  point  numbers  link  the  GRID  output 
points  to  the  difference  quotient  computations. 

4-14.  BACKGROUND.  See  the  discussion  of  difference  quotients  in 
the  VARVARYl  documentation  (Chapter  5). 

4-15.  LIMITATIONS 

a.  The  current  compiled  version  is  limited  to  20  variables. 

b.  The  current  compiled  version  is  limited  to  a  maximum  of  500 
lines  of  difference  quotient  computations  in  each  difference  quo¬ 
tient  block. 

c.  Since  this  routine  utilizes  GRID  output,  the  limitations 
applicable  to  GRID  apply. 

4-16.  RUN  SETUPS 

a.  To  Develop  an  Absolute  ASCI  Program  File 
OMAP.S  name  of  absolute  element. 

IN  03PR0GTEST.DIFFQU0T. 

IN  03PR0GTEST. PARTIAL. 

IN  element  containing  the  driver  PREPR. 

IN  programs  to  be  tested. 

LIB$*FTN8. 

END 
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b.  To  Execute 


@USE  10,  name  of  file  produced  by  GRID  on  unit  10. 

@ASG,A  name  of  above  file. 

@USE  12,  name  of  file  produced  by  GRID  on  unit  12. 

@ASG,A  name  of  this  file. 

@XQT  name  of  absolute  deck  created  earlier. 

[Input  deck] 

c.  Sample  Input  Deck  and  Description 

(1)  Line  1. 

2 

The  number  of  the  variable  for  which  difference  quo¬ 
tients  are  to  be  computed,  in  13  format. 

(2)  Line  2. 

(13F6.0) 

The  format  for  reading  each  line  of  input  data. 

(3)  Line  3. 

5.0  0.5  0.5  0.5 

Initial  values  for  each  variable,  where  variable  1  is 
leftmost. 

(4)  Line  4. 

18.00  9.00  48.00  0.9 

Terminal  values  for  each  variable,  variable  1  leftmost. 

(5)  Line  5. 

3.5  1.5  10.0  0.1 

Step  values  for  each  variable. 
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{1X,8F13.5) 

Format  for  one  line  of  difference  quotient  printouts. 
(7)  Line  7. 

(1X,8(5X,I6)) 

Format  for  printing  out  the  point  numbers  whose  differ¬ 
ence  quotients  are  being  computed. 

d.  Sample  Input  Files  and  Descriptions 

(1)  Unit  10 


1260  4  .0010000 


(6(7X.F13.5)) 

5.00000 

.50000 

.50000 

.50000 

.73821 

5.00000 

.50000 

.50000 

.50000 

.77276 

5.00000 

.50000 

.50000 

.70000 

.89735 

5.00000 

.50000 

.50000 

.80000 

.84199 

5.00000 

.50000 

.50000 

.90000 

.87667 

5.00000 

.50000 

.50000 

1.00000 

.91140 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

This  file  is  in  standard  format;  it  contains  variable  values 
(points)  and  output.  This  file  is  produced  by  GRID  on  unit  10, 
see  GRID  documentation  (Chapter  4)  for  a  more  detailed  descrip¬ 
tion. 


(2)  Unit  12 


5.00000 

.50000 

.50000 

.50000 

.73821 

4.00100 

5.00000 

.50000 

.50000 

.60000 

.77276 

4.00100 

5.00000 

.50000 

.50000 

.70000 

.80735 

4.00100 

5.00000 

.50000 

.50000 

.80000 

.84199 

4.00100 

5.00000 

.50000 

.50000 

.90000 

.87667 

4.00100 

5.00000 

.50000 

.50000 

1.00000 

.91140 

4.00100 

5.00000 

.50000 

10.50000 

.50000 

7.44301 

3.00100 

5.00000 

.50000 

10.50000 

.60000 

7.41817 

4.00100 

5.00000 

.50000 

10.50000 

.70000 

7.39330 

4.00100 

-This  .fU.e  Js.  urocUiced  by  GRID,,  on  unit  J2  and  is  used  by  DIFFQUOT 
to  determine  the  periodicity  of  each' Variable.^'ToV  ff Tiiore' de- *  • 
tailed  description,  see  the  GRID  documentation  (Chapter  7). 


4-31 


CAA-D-80-1 


e.  Sample  RUN  SETUP 


@USE  10,03MAT5. 

@ASG,A  03MAT5. 

@USE  12,03MAT7. 

@ASG,A  03MAT7. 

(aXQT  03PR0GTEST.DIFFQU0T 
1 


(13F6.0) 

5.0  0.5  0.5 

0.5 

18.00  9.00  48.00 

0.9 

3.5  1.5  10.0 

0.1 

(1X,8F13.5) 

(1X,8(5X,16)) 

4-17.  DESCRIPTION  AND  SAMPLE  OUTPUT 

. DIFFERENCE  QUOTIENTS  FOR  VARIABLE  1 

AND  POINTS 


1 

253 

505 

757 

1009 

.00000 

.00032 

.00000 

.00023 

.00015 

.00000 

.00018 

.00011 

.00008 

.00000 

.00015 

.00009 

.00007 

.00005 

.00000 

The  following  figure  (4-1)  may  help  explain  the  output: 

The  slope  of  the  line  to  point 


From 
poi  nt 

1 

253 

505 

757 

1009 

/I 

.0032 

.00023 

.00018 

.00015 

253 

505 

757 

1009 

Slope 

/- 

of  the  line 

.00015 

.00011 

.00008*^ 

.00009 

.00007 

.00005 

.from  1 
to  poi 

3oint-nOi  1-- 
int  no.  253 

Slppe  of  the 
point  no.  25 
no.  757 

‘line  from 
>3  to  point 

Figure  4-1.  Explanatory  Figure 
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For  more  details,  examine  the  description  of  difference  quotients 
in  the  VARVARYl  documentation  (Chapter  2,  Section  II). 

4-18.  DIFFQUOT  ROUTINE  LISTING 


P»R»MCTE:R  IN^ILc  :3,OuTF'TL  =  E  ,r40PTS  =  EIJlltNC  V»LS  =  20.FILEi  =  XO 
PARAMETER  F  IL  £3  =  1 2  •  NC  V  *1- irNO  V  AL^-f  1 ,  fjO  V  AL  2  =  NO  V  AL  1 
C 

DIMENSION  FSTVALfNC  VALTI tLSTVALINOVALSI , STEPS • NO V ALS  I 
DIMENSION  INDEX fNCVALS )»VALUINCVAI.D»»CCORDSINDPTS  fNOVALl ) 
DIMENSION  C IFF (NOV ALS I , FO INT S ( NOP TS ( • CLOCK (NO V ALS  » 
DIMENSION  FERICCtNOVALSJ 
C 

real  EPSL 0N,FSTVAL,L ST VAL, steps  tV ALU* oiff, coords 

c 

INTECER  InSEX.POINTS  •‘lOVAR*"fNDVPMl*NCVSFl*WMICH*VAK,V  ARl 
INTEGER  LOC*PTCT*I*J*V*OLQCK,Pr^ICL  *  PTP  fNtTMflEfi 
C 

CHARACTER*PU  Em  TR  D , '‘‘•TPPT  .  FMTPR  .1  *  EM -^F IL 
C 

REWIND  FILEI 
REWIND  files 

RFAO(  FILEltlftUAO  I  \UM3Cr*NCVARStErSLON 
iUL4D  FORMATt  I5*IStElC1.7  1 

RfAO(FIL£l*lLUlUJ  FMTFIL 
C 

RFAO( infile  .ioniuu  WHICH 
lOCOU  FCRMAT(T3» 

REAS( INFILE  ilUUlU)  ENTRC 
IWIO  FORMAT(  AGO) 

ReAS( INFILE  fPMTRO I  t F ST V AL « V  I  * V: 1 . NC V A?S ) 

REAO( INFiLf .F«TRD I  ( LS T V AL ( V »  t V: l»NO V AR S I 
REA3(  INFIL'"  t^MT^D  I  t  STEFSt  VI  .  V^l.NOVARS) 

READ(  INFILE  .II’OIOI  FMTFRT 
PEAOC INFILE .lOLin 1  FMTER^ 

C 

NaVRPl^NCVARS*-! 

NOVRP2::r!OVFPl  +  l 
NOVRMS:iiOVARS-l 
DO  100  miNOVARS 
CLOCK (1 1-1 
TNOEXd  |:l 
VALU(  ICFSTVALdI 
100  C0NTIf4Ur 

IF  (  WHICH  .EQ.  NOVARS  >  GOTO  ''HO 
DC  2U0  I-WHICH,N0VRH1 
INCEXddlNDEXd  +  l  » 

21JU  CONTINUE 

INDEX (NCVARSI^WHICH 
20U  LOCzp 

CO  1200  IdfNQVARS 
PERIOD(  ’’1-0 
12UU  CONTINUE 
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DO  13UO  1=1. NUMBER 

FEADIFILE3.FMTFTL »  » V ALU ( J  I  . J= 1 . NO VRP2 I 
PTR=INTI VALUCNOVRP?) » 

IF  «  PFRIOOIFTRl  .FT.  D  )  F ER ICO  I PT ” I  =  M A XU  I  I- 1 . U 
IF  »  PFRIOOCll  .NE.  D  1  GOTO  l^SU 
130n  CONTINUE 
STOP 

1150  DO  1375  I=1.N0VARS 
VALU(II=FSTVALm 
i37S  CONTINUE 
GOTO  1000 

'lOU  DC  40t,  VAR  =  N0VAPS.1«-1 
VAR1=VA9 

IF  «  VALUlINOrxCVAR)  1  .LT.  L S T V AL ( I  NCI X t V A P )  )  )  GOTO  7DU 

CL0CK(INDEX( VAP J  J=1 

VALU( INDEX ( VAR) ) =FS T V A L ( INDEX « V AR ) ) 

C  IF  LCC  IS  7FR0.  HAVE  ALREADY  PRINTED  THE  LAST  SDT 
C  OF  DIFFERENCES 

IF  (  LOC  .EG.  n  )  GOTO  400 
WRlTEIOIJTFTL.irOZUI  UHICH 

lunzo  F0RPAH//,lX.2r(lH-l  .•DTFFERENrr  aUOTTENTS  >^00  VASIAPLE 
1  :5.2X.2U(1H-1  ) 

VRITECOUTFTL. 10025) 

10025  FORHATI 21X. • AND  POINTS*) 

WRITE  (OUTFIL.  F)- TP  R2)  «  PC  INTS  •  J  )  .  J  =  1  .LOC) 

WPITE(0UTFIL. 10023) 
irjl)2e  FORMAT)/) 

CO  500  tn.LOC 
OIFFJDrU. 

00  600  J=1.I 

IF  (  J  .EO.  I  )  GOTO  600 

DirF(J)= (COORDS) I .N0V3F I ) -COORDS CJ.NOVRPl) )/ 

1  (COORDS) I. WHICH) -COORDS (J. WHICH) ) 

600  CONTINUE 

WRITE (OUTFIL. FMTP FT )  ) Cl FF ) J ) . J= 1 . I ) 

SUU  CONTINUE 
LOr=U 

400  CONTINUF 

wr:te)outfil. 10030) 

10030  DQRMAT)///) 

STOP 

7  00  VALU) INDEX ) VARI) )  =  VA LU ) IN OEX ) V AR 1 ) ) tS TF P S ) TNO E X ) V  API )  ) 

CLOCK (INDEX ) VAPl ) ) =CL 3CK ) INDE X ) V API  )  )♦! 

1000  CALL  PRFPR( VALU.VALU(NOVPPl)  ) 

L0C=L0C+1 
POINTS) L0C)=1 
DC  1100  I=1.NCVARS 

P0lNTS(L0C)  =  FrlNT3)LCC  )♦ ( CLOCK ) T )-l ) » PER  IOC ) I ) 

1100  continue 

DO  300  I=1.NCVPP1 
COORDS)LOC.I)=VALU):) 

BOU  CONTINUE 

GOTO  JDD  .... 

END 
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4-19.  DIFFQUOT  ROUTINE  FLOWCHART 


I  Read  in  the 
/number  of  points, j 
/of  variables,  and/ 
f  the  zero  / 
approximator  / 


r  liead  in 
the*T11e 
input 
format  i 


/Read  in  the  number/ 
I  of  the  variable  / 
/for  which  differ-/ 

I  ence  quotients  / 
must  be  / 

computed  / 


Read  in  the 
format  for 
the  input 
data 


/  Read  in  the 
Initial ,  terminal / 
and  ‘.cop  / 

values  / 


/  Read  in  the 
I  format  for  one 
/l  ine  of  differ-  y 
I  ence  quotient  / 
_  output  / 


/  Read  in  the 
I  format  for 
printinq  the 
point  numbers 


I  tnltiJiize  the 
/point  descriptor j 
I  indirect  address/ 
r  inq  and  va lue  / 
vectors _ ] 
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CHAPTER  5 

COMMON  SUBROUTINES 


5-1.  INTRODUCTION.  Two  subroutines  are  used  by  both  the 
POINTCOMP  and  GRID  subsystems: 

a.  PARTL 

b.  PREPR 

Detailed  writeups  of  these  subroutines  follow. 


Section  I.  THE  PARTL  SUBROUTINE 

5-2.  INTRODUCTION  TO  THE  PARTL  SUBROUTINE.  This  routine  computes 
the  partial  derivative  at  a  point  numerically  by  using  the  defini¬ 
tion  of  the  right  partial  derivative.  PARTL  is  called  by  POINT¬ 
COMP  and  other  routines  repetitively  to  construct  the  gradient, 
but  it  may  be  used  by  anyone  as  a  standalone  subroutine.  PARTL 
calls  PREPR  which  must  be  a  user-provided  driver  subroutine  whose 
function  is  to  obtain  an  output  value  from  the  program  being 
tested. 


5-3.  BACKGROUND.  Given  a  real  valued  function  f{xj,  x  ), 
the  jth  right  partial  derivative  of  f  at  (xi,  ...,  x„)  is  defined 
to  be:  ^  ” 


1  im„  ' 

h->0+ 


.  x^) 
f(xi. 


Xj-^»Xj->-h»Xj4.^,  •♦♦»Xp)  -  f(x^ . Xp) 

h 


5-4.  DISCUSSION  OF  METHODOLOGY 


a.  The  methodology  is  derived  directly  from  the  definition: 

For  h  =  1/2,  we  approximate  the  jth  partial  derivative  at 
(xj . Xj^)  to  be: 


f(xi,  ...,  Xj-i^Xj+1/2,  Xj^i,  ...,  x^)  -  f(xi . Xq) 

1/2 

b.  We  repeat  this  procedure  for  h  =  1/4,  1/8,  1/16,  etc.  When 
the  difference  between  two  successive  approximations  is  suffi¬ 
ciently  small  (as  determined  by  an  input  parameter),  the  process 
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is  terminated  and  the  approximation  is  returned  as  the  value  of 
the  partial  derivative.  The  process  terminates  automatically 
after  50  approximations. 

5-5.  LIMITATIONS 

a.  As  currently  compiled,  all  inputs  to  PARTL  must  be  real. 

b.  The  function  whose  partials  are  being  computed  should  be 
differentiable. 

5-6.  CALLING  SEQUENCE 

Call  PARTL  (II,  12,  Rl,  R2,  R3,  R4,  R5) 

where 

11  is  an  integer  input  variable  containing  the  number  of  vari¬ 
ables. 

12  is  an  integer  input  variable  containing  the  index  value  of 
the  variable  whose  partial  is  to  be  found.  The  index  value 
refers  to  the  subscript  locating  the  variable  in  the  input 
array  Rl. 

Rl  is  a  real  input  array  containing  the  values  of  the  vari¬ 
ables  at  the  point  at  which  the  partial  is  to  be  evaluated. 

R2  is  a  real  input  variable  containing  a  positive  number.  Any 
number  smaller  than  this  number  will  be  considered  to  be 
zero. 

R3  is  a  real  output  variable  into  which  the  approximate  value 
of  the  partial  will  be  placed  by  PARTL. 

R4  A  real  output  variable  into  which  the  last  increment  tested 
will  be  placed  by  PARTL. 

R5  is  a  real  input  variable  containing  the  value  returned  by 
PREPR  for  the  variable  values  in  Rl. 
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5-7.  PARTL  SUBROUTINE  LISTING 


subroutine  PARTLlNOVARS.WH|CM.V«l..CpS.P«RTVLt INCR.BaSEI 
PARANETCR  N0VALS>20, lNP|LC*S.OUTP|LaA 

dimension  valinovalsi >vali (Novalsi 
INTC6ER  NOVARSiNHiCHtFiRSTtI (J 

RCAL  VAUiVALI  1 1  NCR  i  VALUE  I  .  VALUES  ,PARTYL  iCPS  iBASE 

FIRST-I 
|NCR«|  • 

DO  lOO  laitNOVARS 
VALI  I  I  taVALI  I  I 
loo  CONTINUE 

00  200  Ja| (SO 
|NCRaINCR«l<SI 
VAL|i«HICH|aVALIWHlCHI«iNCN 
CALL  PREPRl VALI .VaLUEII 
IF  IFIRST  (NE.  ll  Goto  SOO 

PART YL« I  Value  1 -base  I /inch 

FIRSTaO 

OOTO  200 

300  IF  UBSH(VALUEl-BAS£l/fNCR>-PAHTYL»  .LT,  EPSI  RETURN 
PARTYL««VALUE|-B*SE|/INCR 
2oO  CONTlNVe 
return 

END 
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Section  II.  THE  PREPR  DRIVER  SUBROUTINE 

5-9.  INTRODUCTION.  This  subroutine  must  be  written  by  the  user 
to  provide  an  interface  between  the  PROGTEST,  PARTL,  VARVARYl, 
GRID,  and  DIFFQUOT  routines  and  the  routine  being  tested.  PREPR 
receives  variable  values  from  POINTCOMP  or  the  other  routines  and 
returns  the  output  value  determined  by  calling  the  program  being 
tested  with  the  given  variable  values  as  input. 

5-10.  DISCUSSION.  The  input  variable  values  are  passed  to  PREPR 
in  array  V,  in  the  same  order  as  the  values  read  in,  i.e.,  the 
leftmost  variable  defined  in  the  input  is  in  V(l),  etc. 

5-11.  PREPR  LAYOUT 

SUBROUTINE  PREPR(V,VALU) 

PARAMETER  NOVALS  =  20 
DIMENSION  V( NOVALS) 

REAL  V.VALU 


Pass  input  values  given  in  V  to  the  program  being  tested. 

Call  program  being  tested  as  a  subroutine  using  input  values 
passed  in  V. 

VALU  =  value  returned  by  the  program. 

END 


i 

i 
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